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APPENDIX  1 


Development  of  Warning  Model 


A  model  to  determine  the  expected  effectiveness  of  a  community 
civil  defense  warning  system  is  described  in  this  appendix. 


I.  THEORETICAL  DEVELOPMENT  OF  MODEL 

The  following  definitions  are  required: 

E  The  event  that  at  least  one  person  in  a  household  or  business 
establishment  is  listening  to  a  radio  or  television  set  with 
enough  attentiveness  so  that  he  would  hear  and  understand 
the  meaning  of  a  warning  if  one  were  transmitted. 

F  The  event  that  no  one  in  a  household  or  business  establish¬ 
ment  is  listening  to  a  radio  or  television  set,  but  there  is 
at  least  one  person  in  the  household  or  business  establish¬ 
ment  who  has  available  a  radio  or  television  set  which  could 
be  used  to  receive  messages. 

N  The  event  that  at  least  one  person  in  a  household  or  business 
establishment  receives  an  alert  from  a  NEAR  installation. 

S  The  event  that  at  least  one  person  in  a  household  or  business 
establishment  is  within  reception  range  of  a  civil  defense 
siren, 

W  The  event  that  a  person  receives  a  warning  message  or  signal, 
believes  it,  and  acts  accordingly. 

Events  E,  N.  F,  and  S  are  assumed  to  be  independently  distributed. 

The  total  population  is  thought  of  as  being  in  one  of  three  locations- 
(a)  in  households,  (b)  in  business  establishments,  or  (c)  m  transit. 
Households  include  all  domiciles.  A  business  establishment  includes 
all  places  of  work  and  also  schools  and  churches.  People  in  transit  are 
those  who  are  driving,  walking,  in  playbrounds,  and  all  those  who  for  one 
reason  or  another  are  not  in  a  household  or  business  establishment. 
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It  is  assumed  that  if  a  person  receives  a  believable  warning  he  will 
spread  the  word  to  every  person  in  his  household  or  business  establish¬ 
ment.  It  is  assumed  that  all  persons  alerted  in  this  manner  would  also 
believe  in  the  validity  of  the  warning  and  act  accordingly. 

The  requirements  that  must  be  met  in  order  to  have  the  event  W 
occur  may  be  specified  in  a  number  of  ways.  Two  different  ways  of 
specifying  these  conditions  are  considered  and  a^e  designated  W ^  and 
W2. 


JLet  W  take  place  if, 

(a)  A  person  hears  a  warning  over  television  or  radio,  or 

(b)  A  person  hears  a  warning  signal  from  a  siren  and  from 
a  NEAR  device,  or 

(c)  A  person  hears  a  warning  signal  from  a  siren  and  then 
turns  on  his  radio  or  television  set  and  receives  a  con¬ 
firmation  of  the  siren  warning,  or 

(d)  A  person  hears  a  warning  signal  from  a  NEAR  set  and 
then  turns  on  his  radio  or  television  set  and  receives 
a  confirnriation  of  the  NEAR  warning  signal. 

These  conditions  may  be  stated  mathematically  as  follows: 

Wj  -  f(E)  U  (SN)  U  (SF)  U  (NF)  J 

It  is  assumed  that  if  a  person  hears  either  a  signal  from  a  siren  or 
a  NEAR  set  and  if  he  has  a  radio  or  television  set  available,  he  will  turn 
on  the  radio  or  television  set  to  ascertain  the  meaning  of  the  siren  or 
WEAR  signal. 

The  probability  of  being  warned  under  these  conditions  is  given  by:^ 

p(Wj)  p(E)  t  (Aj  -  ♦  A^)  (1  -  p(E)) 

where:  Aj  p(NF)  ^  p(NS)  ♦  p(SF) 

A2  p(SN)p(FN)  ♦  p(SF)p(SN)  ♦  p(SF)p(FN) 

Aj  -  p(SN)p{FN)p'SF) 

^  This  expression  is  derived  in  Appendix  la. 
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Let  W  take  place  if, 

(a)  A  person  hears  a  warning  over  television  or  radio,  or 

(b)  A  person  ht^ars  a  warning  signal  frcxn  a  siren,  or 


(c)  A  person  hears  a  warning  signal  from  a  NI^AR  set. 

These  conditions  may  be  stated  mathematically  as  follows; 

W,  -  (E  IJ  S  U  N) 

The  probability  of  being  warned  under  these  conditions  is  given  by; 
p(W2)  =  B,  -  +  Bj 

where ; 

=  p(E)  +p(S)  +p(N) 

B2  =  p(EN)  +  p(ES)  +  p(NS) 

B3  =  p(ENS) 

II.  APPLICATION  OF  THE  MODEL 

The  p^W)  expressed  by  the  model  may  be  interpreted  as  the  expected 
propor  tion  of  the  population  that  will  receive  a  believable  warning.  The 
model  may  be  used  to  assess  the  effectiveness  of  civil  defense  policies 
to  install  warning  sirens  or  NEAR  units  in  a  community  under  various 
assumed  conditions  required  for  a  warning.  The  requirements  for  Wj 
may  be  thought  of  as  assumptions  applicable  when  the  public  expectation 
of  nuclear  war  is  relatively  low  The  requirements  for  Wjn  on  the  other 
hand,  represent  an  ideal  state  of  awareness  which  probably  can  never  be 
reached. 

A  telephone  survey  was  conducted  in  the  City  of  Stamford  with  the 
consent  of  the  Stamford  Office  of  Civil  Defense  to  determine  values  of 
the  parameters  p(E)  and  pfEl  required  in  the  ntodet.  The  telephone 
calls  were  made  between  the  hours  of  ^:i0  -  1 1  :i0  a.  m.  and  1 :30  - 
3:30  p.  m.  ,  on  weekdays  from  the  IJth  \o  the  ^Oth  of  August.  The 
object  was  to  find  typical  values  for  the  working  hours  of  a  summer 
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day  in  Stamford.  On  the  basis  of  468  answered  telephone  queries  in 
which  the  disposition  of  over  2,  000  people  was  directly  or  indirectly 
determined,  estimates  for  p(E)  and  p(F)  were  found  to  have  the  follow¬ 
ing  values: 

p(E)  =  .23 
p(F)  =  .20 

These  figures  may  be  interpreted  as  follows.  Of  the  population  in 
Stamford  not  in  transit,  23%  are  in  a  household  or  business  establish¬ 
ment  where  someone  is  listening  to  a  radio  or  television  set  with  suffi¬ 
cient  attentiveness  to  hear  and. believe  a  disaster  warning  if  one  v;ere 
broadcast.  Of  the  people  in  Stamford  who  are  in  a  household  or  business 
establishment  where  no  one  is  listening  to  a  radio  or  television  set,  20% 
of  them  have  a  television  or  radio  in  the  immediate  vicinity  that  could  be 
turned  on. 

Using  the  estimates  of  p(E)  and  p(F)  established  by  the  survey,  the 
model  was  used  to  calculate  p(Wi)  and  p{W2)  for  various  hypothetical 
values  of  p(S)  and  p(N).  The  results  are  presented  in  Tables  1  and  2. 

The  effectiveness  of  the  warning  system  is  a  strong  function  of  time 
since  the  location  of  people  with  respect  to  warning  devices  varies  con¬ 
siderably  throughout  time.  The  physical  factors  which  govern  the  range 
and  effectiveness  of  individual  devices  also  depend  upon  the  time  of  day, 
season  of  the  year,  etc. 

To  illustrate  the  importance  of  time  as  a  variable,  data  was  taken 
trom  a  study  conducted  in  Chicago^  to  estimate  p(E)  and  p(F)  for  the 
hours  10:30  p.m.  and2K)0a,  m. 

The  calculations  and  assumptions  used  to  estimate  p(E)  and  p(F) 
from  the  data  are  given  in  Appendix  lb.  The  following  values  were 


estimated; 

10 :30  p.  m. 

2  ."00  a.  m. 

P(E) 

.  5?. 

.  06 

P(F) 

.45 

.  87 

proportion  of 

.  07 

.  02 

population  in  transit 

Using  these  values  in  the  model,  p(Wj)  and  p(W2)  were  calculated  and 
are  presented  in  Tables  3,  4,  5,  and  6. 


1  Raymond  W.  Mack,  George  W.  Backer,  The  Occasion  Instant.  The 
Structure  of  Social  Responses  to  Unanticipated  Air  Raid  Warnings. 
Disaster  Study  No.  5,  Publication  945,  National  Academy  of  Sciences- - 
National  Research  Council,  Washington  25,  D.  C.  ,  1961. 
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Table  1 


Probability  of  Being  Warned^  p(Wj),  on  a  Typical  Summer 
Working  Day  in  Stamford,  Connecticut,  as  a  Function  of  the 
Percentage  of  the  Population  which  is  Alerted  by  a  Siren 
and  a  NEAR  Device. 
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100% 
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1.  00 

80% 
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.59 
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.  90 

60% 

.  32 

.43 

.  53 

.  62 

.  70 

.  78 

40% 

.29 

.  37 

.  45 

.  53 

.59 

.66 

20% 

.  26 

.  32 
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0% 
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Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


Table  2 

Probability  of  Being  Warned^  p(W2),  on  a  Typical  Summer 
Working  Day  in  Stamford,  Connecticut,  as  a  Function  of  the 
Percentage  of  the  Population  which  is  Alerted  by  a  Siren 
and  a  NEAR  Device. 
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100% 

1.  00 

1.00 

1.00 

1.  00 

1.  00 

1.00 

80% 

.  85 
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.  91 

.  94 

.  97 

1.00 

60% 

.  69 

.  75 

.  82 

.  88 

.  94 

1.00 

40% 

.  53 

.  63 

.  72 

.  82 

.  91 

1.00 

20% 

48 

.51 

.63 

.  75 

.  88 

1.  00 

0% 

23 

.48 

.  53 

.69 
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1.  00 

0% 

20% 

40% 

60% 
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100% 
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Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


1 


People  in  transit  not  included. 
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Table  3 


Probability  of  Being  Warned^  p(W^i)  10;30  p.  m.  (Chicago) 
as  a  Function  of  the  Percentage  of  the  Population  which  is 
Alerted  by  a  Siren  and  a  NEAR  Device. 


Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


Table  4 

Probability  of  Being  Warned^  p(W|)  at  2K)0  a.  m.  (Chicago) 
as  a  Function  of  the  Percentage  of  the  Populs.tion  which  is 
Alerted  by  a  Siren  and  a  NEAR  Device. 


Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


^  The  7%  of  the  population  estimated  to  be  in  transit  is  not  included. 
^  The  2%  of  the  population  estimated  to  be  in  transit  is  not  included. 
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Table  5 


Probability  of  Being  Warned^  10:30  p.  m.  (Chicago) 

as  a  Function  of  the  Percentage  of  the  Population  which  is 
Alerted  by  a  Siren  and  a  NEAR  Device. 


100% 

1.  00 

1.00 

1.  00 

1.  00 

1.00 

1.00 

80% 
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Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


Table  6 

Probability  of  Being  Warned^  pfWo^  2:00  a.  m.  (Chicago) 
as  a  Function  of  the  Percentage  of  the  Population  which  is 
Alerted  by  a  Siren  and  a  NEAR  Device. 
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Percent  of  Population  Alerted 
by  a  NEAR  Device  p(N) 


^  The  7%  of  the  population  estimated  to  be  in  transit  is  not  included. 
2  The  2%  of  the  population  es*^imated  to  be  in  transit  is  not  included. 


1-7 


Appendix  la 


By  definition: 

Wj  =  [(E)  U  (NF)  U  (NS)  U  (SF)] 
then:  p(W^)  =  p(E)  +  (Aj  -  A2  +  A3)  (1  -  p(E)) 

where:  Aj  =  p(NF)  +  p(NS)  +  p(SF) 

A2  =  p(NS)p(NF)  +  p(SF)p(NS)  +  p(SF)p(NF) 
A3  =  p(SN)p(FN)p(SE) 

Proof : 

p(Wi)  =  p  [(E)  U  (NF)  U  (NS)  U  (SF)] 
call  [(NF)  U  (NS)  U  (SF)]  =  B 
then  by  substitution: 

p(Wj)  =  p(E  U  3) 

"  p(E)  +  p(B)  -  p(E)p(B) 

=  P(E)  +p(B)  [1  -  p(E)] 

It  remains  to  be  shown  that: 

P(B)  =  Aj  -  +  A3 

By  definition: 

B  =  [(NF)  U  (NS)  U  (SF)J 

then:  p(B)  =  p  [(NF)  U  (NS)  U  (SF)  j 

call  [(NS)  U  (SF)]  =  C 

then:  p(B)  ^  p  [(NF)  U  (C)  J 

=  p(NF)  +p(C)  -  p(NF)p(C) 
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But:  p(C)  =  P(SN)  +  p(SF)  -  p(SN)p(SF) 

Substituting  this  expression  for  p(C)  in  above: 


p(B)  =  p(NF)  +  [p(SN)  +p(SF)  -  p(SN)p{SF)]  - 

p(NF)  [p(NS)  -f  p(SF)  -  p(NS)p(SF)] 
Rearranging  terms: 

p(B)  =  p(NF)  +  p(NS)  +  p(SF)  -  p{NS:.p(NF)  - 

p(SF)p(NS)  -  p{SF)p(NF)  +  p(SN]p(FN)p{SE) 

=  Aj  -  A2  + 


1-9 


Appendix  lb 


The  following  values  were  taken  from  The  Occasion  Instant, 
pages  29  and  30: 

What  people  were  doing  at  10:30  p.  m. 

42%  Watching  television 
13%  Conversing 
13%  Preparing  for  bed 
8%  Listening  to  the  radio  or  reading 


Remaining : 


Household  chores 
Driving 
Walking,  etc. 


All  of  the  activities 
less  than  seven  percent  ea 


100% 

From  the  above  values,  the  following  more  detailed  breakdown  was 
constructed: 

42%  Watching  television 
13%  Conversing 
13%  Preparing  for  bed 
4%  Reading 
4%  Listening  to  the  radio 

5%  Household  chores,  sleeping,  miscellaneous 

81% 

7%  Driving  and  walking 
12%  Working 

100% 

p(E)  was  calculated  as  follows; 

p(E)  =  .42  Proportion  watching  television 

+  ,04  Proportion  listening  to  radio 

+  .  06  Half  of  the  working  population 

=  .  52  at  10  ;30  p.  m. 
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p(F)  was  calculated  as  follows: 

From  I960  census  data: 

,  915  proportion  of  households  have  a  radio 
.  863  proportion  of  households  have  a  television 

Assuming  independent  distribution  of  radios  and  televisions,  the 
probability  that  a  household  will  have  a  radio  or  a  television  set  is 
given  by : 


each 


p(RUT)  =  .915  +.  863  -  915){.  B63) 

=  .98 

This  means  approximately  98%  of  the  households  do  have  either  a 
radio  or  television  set. 


From  above,  42%  plus  4%  are  listening  to  radio  or  television. 
This  leaves  54%  who  are  not  listening  to  radio  or  television.  Out 
of  this  54%,  7%  are  in  transit,  leaving  47%.  p{F)  is  now  arbitrarily 
taken  to  be  a  high  percentage  of  this  remainder. 


p(F)  =  ,  45  at  10:30  p.m. 

Parameters  at  2:00  a.  m. 

Estimate  of  what  people  are  doing  from  Chicago  data: 

12%  Working 

2%  Driving  or  walking  (in  transit) 

86%  At  home 

p(E)  =  No  listeners  at  home  and  half  of  the  workers 
=  0  +  .  06 
p(E)  =  .  06 

p(F)  =  98%  at  home  have  radio  or  television.  Say  half  of 

working  population  not  listening  have  radio  or  tele¬ 
vision  available 

p(F)  =  (.  98)(.  86)  f  1/2  .06 

p(F)  =  .  87 
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Development  of  Assignment  Model 


I.  INTRODUCTION 

In  any  community  it  is  likely  that  there  will  be  some  areas  where  there 
are  more  people  than  shelters,  and  other  areas  where  there  are  more  shelters 
than  people.  It  is  also  likely  that,  unless  something  is  done  to  prevent  it, 
excessive  and  unnecessary  casualties  will  result  from  overcrowded  conditions 
in  some  shelters  while  others  remain  v^acant  and  from  the  exposure  of  people 
to  lethal  doses  of  radiation  while  searching  aimlessly  for  available  shelter. 

In  this  appendix  a  technique  is  presented  for  assigning  persons  from  areas 
of  excess  population  to  areas  of  excess  shelter  in  such  a  way  as  to  minimize 
the  over-all  distance  traveled  and  to  equalize  the  degree  of  overloading,  if 
any,  on  the  shelters  in  the  community.  This  is  done  on  the  assumption  that 
distance  is  proportional  to  time.  (It  is  desired  to  minimize  travel  time  but 
not  practical  to  do  so  with  available  data  at  this  time. )  It  is  also  assumed 
that  under  actual  emergency  conditions  people  will  not  be  turned  away  from 
any  shelter  which  can  possibly  accommodate  them. 

Distance  was  selected  as  a  criterion  rather  than  time  to  shelter  primarily 
because  of  the  greater  confidence  which  could  be  placed  in  its  measurement. 
The  existing  network  of  roads  and  topography  of  the  area  make  a  limited  num¬ 
ber  of  routes  from  origin  to  destination  feasible,  and  the  distance  which  must 
be  traveled  can  therefore  be  estimated  with  precision  and  confidence.  The 
time  which  would  be  required,  on  the  other  hand,  would  be  affected  by  whether 
people  walked,  drove,  or  were  taken  to  shelter,  and  by  the  type  and  extent  of 
traffic  control  and  direction- -factors  which  are  highly  variable  and  subject  to 
change  by  decree.  It  was  felt  best,  therefore,  to  assign  people  on  the  basis 
of  distance,  and  determine  the  times  required  to  traverse  that  distance  as  a 
function  of  transportation  policies  in  a  later  analysis.  The  means  of  deter¬ 
mining  the  time  required  to  reach  shelter  under  different  transportation 
conditions  is  discussed  in  Appendix  3. 

H.  PROCEDURE 


The  procedure  incluciis;  (1'/  the  collection  of  data  on  the  population  of 
the  community,  the  number  of  shelters  in  it,  and  the  various  routes  which 
might  be  used  by  people  seeking  shelter  :  (2)  the  calculation  of  an  overload 
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factor  where  the  total  population  exceeds  the  total  number  of  shelters,  the 
adjustment  of  the  number  of  shelters  in  each  tract  using  this  overload  factor, 
and  the  determination  of  the  number  of  excess  people  or  excess  shelters  for 
each  area;  and  (3)  the  assignment  of  persons  from  areas  of  excess  popula¬ 
tion  to  areas  of  excess  shelter  capacity,  using  rules  which  assure  at  least 
a  "good"  solution  to  the  problem.  These  three  steps  are  discussed  more 
fully  under  separate  headings  below. 

A.  Data  Sources 


The  principal  sources  of  data  used  in  this  analysis  were  the  I960 
Census  of  the  United  States  (6),  the  final  report  of  Phase  II  of  the  National 
Fallout  Shelter  Survey  (2),  and  a  map  of  the  City  of  Stamford,  Connecticut. 
Because  the  Census  Bureau  data  were  reported  by  census  tracts  and  the 
shelter  survey  data  were  reported  by  National  Location  Code  locations,  it 
was  necessary  to  get  the  correspondence  between  the  two.  This  information 
is  available  for  all  regions  in  a  series  of  volumes  prepared  for  OCD  and 
OEP  by  the  Bureau  of  the  Census  (5).  The  National  Fallout  Shelter  Survey 
Report  Formats  and  Descriptions  of  Content  (3)  was  found  to  be  indispensable 
in  interpreting  the  data  in  the  report  itself. 

Population.  Assignments  were  based  on  resident  population  but 
could  equally  well  have  been  based  on  daytime  or  nighttime  population. 
Resident  populations  are  given  in  references  2,  5,  and  6,  but  those  in 
reference  6  are  suspect,  a  machine  error  having  been  found  in  some  of 
the  readouts  (4).  Daytime  and  nighttime  populations  are  given  only  in  the 
survey  report  (2)  and  then  only  for  those  locations  which  have  some  poten¬ 
tial  shelter  spaces.  These  can  be  estimated,  where  necessary,  from  other 
data  contained  in  the  census  reports  (6). 

Shelters.  In  the  Stamford  assignments,  total  possible  shelters 
with  a  protection  factor  of  100  or  greater  (after  improvements)  were  used. 
The  data  were  taken  from  reference  2.  where  they  were  found  in  "Summary 
A"  for  each  individual  tract  under  the  heading  "TOT  POSS  SPACES.  "  (See 
reference  3,  page  S,  row  24  under  column  23.)  Depending  on  the  purpose 
to  be  served,  data  for  total  existing  spaced  w’lth  a  protection  factor  of  100 
or  better  or  for  total  existing  spaces  with  a  protection  factor  of  40  or  better 
could  have  been  used. 

Routes  to  Shelter.  Maps  of  the  city  were  obtained  from  the  city 
engineering  department.  These  were  readily  available  for  all  cities  con¬ 
tacted  and  provided  sufficient  detail  for  the  analyst  to  draw  in  the  tract 
boundaries,  pick  out  by  eye  a  reasonable  route  from  cai  h  origin  to  each 


destination,  and  measure  the  distance  along  each  such  route.  There  will  be 
times,  of  course,  when  unique  local  conditions  which  are  not  clear  on  the 
map  will  make  a  choice  a  bad  one,  and  the  services  of  a  local  person  could 
be  used  to  advantage.  This  situation  was  not  encountered  in  the  selection 
of  shelter-seeking  routes  for  Stamford,  Connecticut. 

B.  Data  Processing 

Overload  Factor.  In  many  communities  there  will  be  fewer  total 
shelters  than  persons.  For  these  communities  an  overload  factor  is  com¬ 
puted.  It  is  the  ratio  of  the  total  population  to  the  total  number  of  shelters, 
or,  the  average  number  of  persons  per  shelter.  This  factor  is  applied  to 
the  number  of  shelters  in  each  tract  to  get  a  new,  "adjusted,  "  number  of 
shelter  spaces.  For  communities  with  at  least  as  many  shelter  spaces  as 
people,  the  adjustment  factor  is  1.0.  All  assignments  are  based  on  these 
adjusted  numbers  of  shelter  spaces. 

Shelter /Population  Excesses.  For  each  tract  the  difference  between 
the  adjusted  number  of  shelter  spaces  and  the  tract  population  is  taken  to  get 
the  number  of  excess  persons  or  the  number  of  excess  shelters  for  that  tract. 
Tracts  with  excess  population  then  become  origins,  and  tracts  with  excess 
shelters  become  destinations. 

III.  SETTING  UP  THE  PROBLEM 


Table  1  presents  the  basic  data  for  the  problem  in  matrix  form.  There 
are  sixteen  areas  of  excess  population  and  eight  areas  of  excess  shelters, 
for  a  total  of  twenty-four  areas,  m  this  particular  case,  the  number  of 
excess  shelters,  44,  397,  is  exactly  equal  to  the  number  of  excess  persona. 
Were  this  not  so,  a  hypothetical  seventeenth  origin  would  be  added  to  the 
matrix  to  make  the  number  of  persons  equal  to  the  number  of  shelters. 

The  excess  populations  and  shelter  spi  ~es,  where  they  exist,  are 
shown  in  the  marginal  row  and  column  of  Table  1.  The  distance  in  arbi¬ 
trary  units  from  each  area  of  excess  population  to  each  area  of  excess 
shelter  capacity  is  given  in  the  body  of  the  table.  The  problem  is  to  assign 
the  excess  population  to  the  excess  shelters  in  such  a  way  as  to  minimise 
the  total  number  of  man-units  traveled  in  reaching  shelter.  This  is  the 
standard  "Transportation  Problem  which  is  of  the  form;^ 


^  See  Gass  (reference  1,  pages  6-81  or  any  standard  text  on  linear 
programming  for  a  complete  statement  of  this  problem. 
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Table  1 


Stamford,  Connecticut 
Distance  From  Areas  of  Excess  Population 
to  Areas  of  Excess  Shelter 
(Arbitrary  Units) 


Origin 

1 

11 

12 

.  Destination 

16  17 

18 

22 

23 

Excess 

Population 

2 

143 

140 

122 

136 

150 

146 

163 

163 

1,396 

3 

157 

139 

120 

150 

164 

160 

177 

183 

2,646 

4 

104 

85 

73 

97 

111 

107 

124 

134 

2,464 

5 

83 

63 

31 

56 

70 

66 

83 

67 

2,864 

6 

64 

44 

31 

57 

71 

67 

c^4 

87 

3,768 

7 

118 

61 

58 

107 

117 

103 

132 

121 

2,355 

8 

36 

28 

13 

52 

62 

48 

77 

85 

2,559 

9 

73 

32 

39 

62 

72 

58 

87 

104 

1,958 

10 

58 

21 

26 

47 

57 

43 

72 

86 

2,  180 

13 

22 

29 

12 

15 

29 

25 

42 

42 

2,  987 

14 

21 

55 

47 

32 

43 

36 

23 

10 

1,234 

15 

6 

39 

41 

20 

31 

21 

11 

9 

4,  351 

19 

35 

41 

60 

50 

31 

22 

37 

53 

4,535 

20 

41 

36 

56 

47 

26 

17 

37 

54 

2,  813 

21 

16 

38 

44 

34 

9 

18 

17 

31 

3,  824 

24 

40 

71 

78 

68 

43 

52 

42 

57 

2,463 

Excess 

Shelters 

29,091 

318 

3,392 

1,914 

2,431 

673 

1,933 

4,645 

44,  39? 

Readings  proportional  to  distance,  derived  from  a  mechanical  device 
used  to  measure  miles  on  -  map. 
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a.  =  excess  population  in  area  i 

bj  =  excess  shelters  in  area  j 

distauice  from  area  i  to  area  j 

number  of  excess  persons  from 
area  i  assigned  to  area  j 

Where  it  is  desired  to  select  Xjj  such  that  L 
to  the  constraints: 


I  a. 

=  L  b. 

1 

j  ^ 

L  X.. 

=  b. 

i  u 

J 

I  X.. 

=  a. 

j 

1 

IV.  SOLUTION  OF  THE  PROBLEM 


X.j  Ci;  is  a  minimum  subject 


Two  totally  different  methods  of  solving  the  problem  were  found  feasible, 
with  their  desirability  depending  on  the  time  available  and  the  degree  of 
accuracy  required.  For  preliminary  checks  and  estimating  purposes,  a 
heuristic  "eyeball"  technique  provided  a  solution  which,  in  this  particular 
case,  was  only  5%  or  so  less  than  optimum,  and  which  took  less  than  one- 
tenth  the  time  required  to  reach  an  optimum  solution.  Such  good  results, 
however,  are  not  to  be  expected  every  time,  and  one  must  not  lose  sight 
of  the  fact  that  this  technique  provides  no  objecrive  means  of  determining 
just  how  good  a  solution  is.  Therefore,  for  more  precise  work,  as  in 
planning  and  making  actual  assignments,  the  time  required  to  reach  an 
optimum  by  more  formal  means  is  not  only  well  spent,  but  the  only  sound 
policy.  For  the  optimum  solution,  a  technique  described  by  Gass  (refer¬ 
ence  1,  pages  145-152)  was  employed.  Although  it  took  longer,  it  was 
simple  enough  that  clerical  help  could  use  it  without  difficulty  to  solve 
the  problem,  and  it  had  the  advantage  common  to  formal  mathematical 
techniques  that  it  measured  its  own  accuracy.  Each  of  these  techniques 
is  discussed  briefly  below. 
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Heuristic  Method.  To  get  an  "eyeball"  solution,  one  scans  the  rows 
and  columns  of  the  matrix  looking  for  a  cell  of  least  cost.  When  such  a 
cell  is  found,  as  many  excess  persons  as  possible  are  assigned  to  that 
shelter  area.  When  no  cells  of  least  cost  remain  to  which  assignments 
can  be  made,  the  matrix  is  scanned  to  determine  the  least  costly  place 
to  assign  the  remaining  excess  persons.  When  all  of  these  excess  per¬ 
sons  have  been  assigned,  the  matrix  is  scanned  to  determine  whether 
any  adjustments  can  be  made  in  the  assignments  to  reduce  the  over-all 
cost.  No  attempt  is  made  to  make  all  such  adjustments,  because  to  do 
so  would  be  to  seek  an  optimum  solution  which  can  much  more  easily 
and  reliably  be  done  using  a  formal  approach:  instead,  only  the  "obvious" 
changes  are  made,  and  a  "good"  rather  than  the  "best"  solution  is  accepted. 

Optimum  Method.  Using  this  method,  one  starts  with  what  is  called 
a  "basic  feasible  solution.  "  This  is  a  solution  which  assigns  all  of  the 
excess  people  to  all  of  the  excess  shelters  in  exactly  m+n-  1  assignments, 
where  m  is  the  number  of  areas  of  excess  shelters  and  n  is  the  number  of 
areas  of  excess  population.  This  need  not  be  a  good  solution  and  can  be 
obtained  by  a  very  simple  technique  called  the  "Northwest  Corner  Rule" 
(reference  1,  page  146).  Starting  with  this  solution,  each  iteration  re¬ 
quires  first  mn  algebraic  sums  of  the  form  a  +b  -  c  to  determine  which 
cells  should  be  adjusted;  then  one  or  two  adjustments  involving  the 
removal  of  persons  previously  assigned  to  one  area  and  their  reassign¬ 
ment  elsewhere  subject  to  the  original  constraints.  For  a  city  the  size 
of  Stamford,  approximately  twenty  iterations  would  be  required,  and  the 
entire  manual  process  would  take  on  the  order  of  two  and  one-half  work¬ 
ing  days.  Computer  programs  are  available  for  this  problem,  and  com¬ 
mercial  computer  organizations  exist  which  will  deliver  a  solution  for 
two  to  three  hundred  dollars  per  matrix. 

The  Solution.  Table  Z  presents  the  solution  to  the  Stamford  problem 
obtained  using  the  optimum  technique:  no  eyeball  solution  is  presented. 

This  solution  is  not  necessarily  unique,  but  it  is  the  best  obtainable.  The 
over-all  number  of  man-miles  traveled  by  those  seeking  shelter  outside 
of  their  own  area  is  101,  109,  and  the  average  individual  distance  traveled 
is  2.  28  miles.  This  means  that  at  3  miles  per  hour  walking  speed  50%  of 
the  excess  persons  will  be  sheltered  within  three-quarters  of  an  hour. 
Assuming  that  all  of  those  who  are  assigned  to  their  own  areas  find  shelter 
within  this  time,  with  no  other  form  of  transportation  involved,  over  75% 
of  the  total  population  of  Stamford  will  be  sheltered  within  45  minutes. 
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Table  2 


Stamford,  Connecticut 

Assignment  of  Persons  From  Areas  of  Excess  Population 
to  Areas  of  Excess  Shelter 


7 

8  1204 

9  1958 

10  266 

13  2987 

14 

15  940 

19  4535 

20 

21  1600 

24  2463 


318  2037 

1355 

1914 
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APPENDIX  3 


Assignment  of  Shelters  by  Census  Tract 
I.  INTRODUCTION 

Within  a  given  civil  defense  area  persons  can  be  rationally  assigned  from 
census  tracts  of  excess  population  to  census  tracts  of  excess  shelter  capacity 
using  the  techniques  of  linear  programming  (2).  Given  these  assignments,  it 
is  possible  to  determine  the  pattern  of  traffic  for  the  community  under  certain 
assumptions  and  to  express  this  pattern  In  terms  of  number  of  persons  in  the 
shelter  area  versus  time.  This  appendix  presents  a  technique  for  accomplish¬ 
ing  this  latter  task,  using  Stamford,  Connecticut,  as  a  case  example. 

The  assumptions  used  in  this  analysis  are; 

Aq!  All  persons  will  travel  by  private  automobile  or  walk; 

Aj:  All  persons  moving  entirely  within  a  census  tract  will 

travel  on  foot; 

A2'  Persons  moving  from  one  tract  t )  another  will  travel  by 
automobile  if: 

a.  The  straight-line  distance  from  the  center  of  the  origin 
tract  to  the  center  of  the  destination  tract  is  greater 
than  one  mile,  and 

b.  there  is  at  least  one  car  in  the  household  unit.  ^ 

A^:  No  more  than  one  car  from  each  household  unit  will  be  used; 

A^:  Should  the  number  of  persons  assigned  to  tracts  beyond  one  mile 

exceed  the  number  of  persons  in  household  units  having  one  or 
more  cars,  such  excess  persons  are  assumed  to  ride  with 
persons  having  automobiles; 

A^:  Traffic  moves  from  origin  to  destination  exactly  as  assigned  by 

the  linear  program; 


^A  "household  unit"  generally  corresponds  to  a  family  unit.  See  reference  3 
for  a  precise  definition  of  tliis  term. 


A^:  All  traffic  from  a  given  origin  •will  travel  in  a  wave  of  constant 

length  determined  by  the  constant  speed  of  the  automobile,  the 
average  length  of  cars  in  the  wave,  the  number  of  cars  in  the 
wave,  and  the  constant  separation  between  cars  in  the  wave; 

Ay:  At  each  intersection,  the  first  wave  to  arrive  will  occupy  the 

intersection  until  it  has  completely  passed.  Other  waves  wait¬ 
ing  at  the  intersection  may  enter  only  when  the  wave  has  passed 
or  v/hen  a  gap  appears  in  it.  When  a  waiting  wave  so  enters  an 
intersection,  it  vill  continue  to  occupy  it  until  it  is  completely 
through  or  until  another  waiting  wave  enters  at  a  gap. 

Ag:  When  part  of  a  wave  must  be  diverted  to  a  proximate  origin, 

that  part  shall  be  the  last  part  of  the  wave. 

A^:  The  effects  of  starting  and  stopping  times  and  road  obstructions 

other  than  moving  cars  will  be  neglected. 

II.  SETTING  UP  THE  PROBLEM:  BASIC  DATA 


Three  basic  items  are  required  to  solve  this  problem:  (1)  a  set  of  assign 
ments  of  excess  persons  who  will  walk  and  of  excess  persons  who  will  drive 
from  origin  tracts  to  destination  tracts;  (2)  information  concerning  the  avail- 
abilit'y  of  automobiles  in  the  origin  tracts;  and  (3)  information  about  the  trans¬ 
portation  routes,  including  the  distances  between  points  on  the  automobile 
routes,  and  the  maximum  and  minimum  distances  walked  by  those  on  foot. 

The  assignment  of  excess  persons  to  tracts  having  excess  shelter  capa¬ 
city  is  accomplished  by  means  of  the  computational  procedure  dis¬ 
cussed  in  Appendix  2.  The  determination  of  who  will  walk  and  who 

will  ride  is  accomplished  by  measuring  the  straight-line  distance  from  center 
of  origin  tract  to  center  of  destination  tract.  Where  the  distance  is  one  mile 
or  greater,  automobiles  are  used;  where  the  distance  is  less  than  one  mile, 
the  excess  population  moves  on  foot. 

The  number  of  cars  departing  from  each  origin  is  determined  by  taking 
from  Reference  3  the  total  number  of  housing  units  having  one  or  more  auto¬ 
mobiles  available  and  comparing  it  with  the  number  of  housing  units  whose 
occupants  will  leave  the  origin.  The  smaller  of  these  numbers  represents 
the  number  of  cars  leaving  the  origin.  That  is,  for  each  origin: 
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Where:  A  is  the  number  of  cars  on  the  road; 

r 

is  the  number  of  occupied  housing  units 
with  automobiles; 

Hq  is  the  number  of  occupied  housing  units; 
Pg  is  the  excess  population; 

Pj.  is  the  resident  popvilation. 


The  number  of  cars  going  to  each  of  several  destinations  from  a  given 
origin  is  assumed  proportional  to  the  number  of  people  assigned  to  each 
destination  from  that  origin.  That  is: 


A.. 

ij 


AjPii 


Where: 


is  the  number  of  automobiles  traveling  from 
origin  i  to  destination  j 


A^  is  the  total  number  of  cars  leaving  origin  i 

Pg  is  the  number  of  excess  people  from  origin  i 
assigned  to  destination  j 

P^  is  the  excess  population  of  origin  i 

As  soon  as  the  assignments  of  excess  persons  to  shelter  areas  is  known, 
transportation  routes  can  be  determined  and  distances  measured.  Walking 
routes  are  measured  from  a  point  on  the  periphery  of  the  destination  tract 
along  the  shortest  possible  path  to  the  nearest  and  then  to  the  most  distant 
point  in  the  origin  tract,  giving  two  figures  which  can  later  be  converted,  at 
any  assumed  walking  rate,  to  times  of  first  and  last  arrival  of  foot  traffic 
from  the  given  origin. 

Automobile  routes  must  be  more  carefully  thought  out.  Traffic  is 
generally  assumed  to  flow  from  the  small  roads  and  lanes  onto  the  larger 
routes  and  main  arteries;  but  exceptions  are  made  where  the  number  of 
intersections  can  be  reduced,  where  a  parallel  flow  of  traffic  can  be 
established,  or  in  general  where  the  rapid  flow  of  traffic  will  be  aided. 
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It  is  necessary  to  measure  the  distance  between  each  two  adjacent  points 
along  the  routes,  where  a  point  is  a  destination,  an  origin,  or  an  intersec¬ 
tion  of  two  routes.  The  over -all  distances  measured  are  between  the  center 
of  the  origin  tracts  and  the  center  of  the  corresponding  destination  tracts. 

III.  SOLUTION  OF  THE  PROBLEM 

The  problem  is  first  solved  for  automobile  traffic  alone,  then  for  foot 
traffic.  These  are  then  combined  to  provide  the  over-all  arrival  rates  and 
times  for  each  destination  tract. 

In  presenting  the  technique  used  in  solving  the  problem,  Stamford, 
Connecticut,  is  used  as  an  example.  The  Assignment  Matrix  for  Stamford 
appears  in  Table  1.  The  necessary  automobile  data  are  shown  in  Table  2. 
The  over-all  flow  of  traffic  is  presented  in  schematic  form  in  Figure  1,  All 
calculations  are  for  the  single  case  of  20  mph  driving  speed,  3  mph  walking 
speed,  and  single  lane  vehicular  traffic. 

A.  Automobile  Traffic 


The  first  step  is  to  determine  the  "wavelength"  of  the  cars  leaving 
each  origin.  This  is  the  interval  of  time,  in  hours,  the  moving  wave  of  cars 
requires  to  pass  a  stationary  point.  It  is  given  by  the  expression: 


Where: 


ns 

5280s 

Wgjj  is  the  wavelength,  in  hours,  at  speed  s, 
of  the  wave  travelling  from  origin  i  to 
destination  j. 


N-  is  the  number  of  cars  in  that  wave, 

L  is  tlic  unit  length  of  the  cars,  in  feet, 

s  is  the  speed  at  which  the  cars  are  travel¬ 

ling,  in  mph. 


n  is  the  number  of  car-lengths  separation 

between  cars  in  the  wave,  in  feet  per 
10  mph  of  speed. 

For  L  16,  S  20,  and  n  1.  6,  this  is  numberically  equal  to  (6.06 
10*"^)  (Njj),  and  Table  3  is  obtained  by  a  series  of  constant  multiplications. 


Table  1 


Stamford,  Connecticut 
Assignment  of  Persons  from  Areas  of 
Excess  Population  to  Areas  of  Excess  Shelter 


2 

3 

4 

5 

6 

1396 

2646 

2464 

2864 

3768 

7 

318 

2037 

8 

1204 

(1355) 

9 

1958 

10 

266 

1914 

13 

(2987)* 

14 

1234 

15 

(940) 

(3411) 

19 

4535 

20 

2140 

(673) 

21 

(1600) 

(291) 

(1933) 

24 

2463 

♦Numbers  in  parentheses  indicate  persons  who  will  walk  because  of  proximity 
of  origin  and  destination  tracts. 
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only  those  persons  assigned  to  shelter  areas  more  than  one  mile  distant. 


Table  3 


Stamford,  Connecticut 
Wavelengths  by  Origin  and  Destination 


s 

=  20  mph;  L  =  16  ft 

. ;  n  =  1,5 

Origin 

Oest. 

No.  of  Cars 

Wavelength 

Rate 

(i) 

(j) 

(Cij) 

(^20ij) 

(persons  per  hour) 

7 

1 

11 

82 

.  050 

6,  385.50 

7 

12 

527 

.  319 

6,  385.50 

8 

1 

298 

.  181 

6,  666. 01 

9 

1 

547 

.  331 

5,907.00 

10 

1 

75 

.046 

5,857.50 

10 

16 

548 

.  332 

5,857.50 

14 

23 

308 

.  187 

6,  616.50 

19 

1 

1147 

.  695 

6,  517.50 

20 

17 

638 

.  387 

5,527.50 

24 

1 

612 

.371 

6,663.00 

1 


Rate 


W.. 

ij 


r 
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ORIGIN 


OCSTtNATION 


Note:  In  such  a  term  as  1.43(.  072).  the  first  figure 
that  is  1.43,  represents  distance  in  miles, 
the  figure  is  parenthesis,  that  is  (.072)  repre 
sents  time  in  hours. 


Figure  1.  Stamford,  Connecticut.  Schematic  Diagram  of 
Automobile  Traffic  Flow.  (Not  to  scale.  ) 
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The  second  step  is  to  convert  all  distances  to  times.  This  has  already 
been  done  in  Figure  1,  where  the  times  appear  in  parentheses  following  the 
distances. 

When  these  two  preliminary  steps  have  been  completed,  the  traffic  flow 
can  be  analyzed  by  inspection,  using  a  worksheet  to  indicate  the  time  of 
arrival  of  the  first  and  last  cars  from  each  origin  at  each  destination. 

The  flow  from  origin  14  to  destination  23,  for  example,  is  simple.  The 
first  car  will  arrive  at  23  at  .  080  hours,  and  the  last  will  arrive  .  187  hours 
later,  or  at  .  267  hours. 

The  flow  from  origins  2,  3,  4,  5,  and  6  is  considerable,  but  is  very 
simple  to  compute.  The  first  car  will  arrive  from  the  nearest  origin  -  -  5  -  - 
at  .  153  hours.  Inspection  of  Figure  1  will  show  that  the  flow  from  each  of 
these  five  origins  will  encounter  one  or  more  delays  on  the  way  to  the  destina¬ 
tion  and  that,  therefore,  there  will  be  no  gaps  in  the  composite  wave  entering 
destination  1.  Thus,  the  arrival  time  of  the  first  car  in  the  composite  wave 
will  be  .  153  hours,  the  interval  between  first  cars  from  consecutive  waves 
will  be  equal  to  the  lengths  of  the  waves,  and  the  time  of  arrival  of  the  last 
car  from  the  last  wave  will  be  .153  hours  plus  the  sum  of  the  five  wavelengths. 

The  flow  from  origins  7,  8,  9,  and  10,  however,  is  both  considerable 
and  harder  to  compute.  Two  complicating  factors  not  previously  encountered 
are  present  in  that  pattern;  first,  the  diversion  of  portions  of  the  composite 
wave  to  proximate  destinations;  and  second,  the  existence  of  an  intersection 
where  traffic  crosses  without  merging.  The  diversions  leave  gaps  in  the 
wave  which  might  be  reduced  or  eliminated  as  the  wave  encounters  further 
delays,  or  which  might  be  carried  intact  all  the  way  to  the  ultimate  destina¬ 
tion.  Where  these  exist,  they  present  the  opportunity  for  waves  already 
waiting  at  an  intersection  to  preempt  the  main  artery  and  alter  the  sequence 
of  arrivals.  The  non-merging  intersection,  on  the  other  hand,  must  be 
examined  to  determine  whether  the  two  waves  using  it  will  attempt  to  do 
BO  at  the  same  time.  If  so,  an  existing  gap  in  one  or  both  waves  could  be 
enlarged  or  compressed;  or,  at  the  very  least,  one  of  the  waves  would  be 
delayed  in  reaching  its  destination. 

Because  this  pattern  illustrates  many  of  the  situations  encountered  in 
this  type  of  analysis,  it  will  be  discussed  in  some  detail. 

Inspection  of  Figure  1  show-s  that  the  first  car  to  reach  destination  1  will 
be  from  origin  10.  It  .’/ill  arrive  at  .  146  (hours)  and  the  75th  will  arrive  at 
.  192.  The  remaining  548  cars  are  diverted  to  16  where  the  first  arrives  at 
.  1  7o  and  the  last  at  .  502. 
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Further  inspection  reveals  that  waves  9  and  7  will  follow,  in  that  order, 
right  on  the  tail  of  wave  10,  with  the  first  82  cars  from  7  going  to  11,  and 
the  last  527  going  to  12.  The  first  car  from  9  will  arrive  at  1  at  .  146  +  .  378 
(the  entire  wavelength  of  10)  =  .524  and  the  last  at  .855  because  the  gap  left 
in  wave  10  by  the  diversion  of  part  of  it  was  unchanged  by  further  delays  to 
wave  10. 

The  last  car  from  9  passed  the  road  to  1 1  at  .  013  +  .  01 1  +  .  005  +  .  022 
+  .  378  +  .  331  =  .  760  and,  therefore,  the  first  would  reach  1 1  at  .  760  +  .  022 
=  782  hours,  and  the  last  at  .  782  +  .050  =  .832  hours.  The  83rd  car  from 

7  passed  the  road  to  12  at  .013  +  .011  +  .  005  +  .  378  +  .  331  +  .  050  =  ,  788  hours, 
and  would  arrive  at  the  non-merging  intersection  .  013  hours  later,  or  at  .801 
hours;  and  assuming  no  delay,  would  arrive  at  12  at  .026  hours  later  at  .82  7 
hours.  The  last  car  from  7  would  then  arrive  at  12  at  1.  148  hours. 

Now,  to  check  the  above  assumption;  the  last  car  from  9  would  pass  the 
point  where  cars  from  8  waited  at  .  765.  Therefore,  the  cars  from  8,  which 
had  been  waiting  since  .  065  to  enter  the  main  road,  would  enter  at  .  765  and 
would  clear  the  non-merging  intersection  at  .  765  +  .  181  -  .026  =  .920  hours. 
From  this  it  follows  that  the  cars  from  7  would  have  to  wait  at  the  non- merging 
intersection  ,117  hours,  so  that  arrivals  at  12  would  be  at  .946  and  1,265 
instead  of  .  829  and  1.  148. 

Since  the  wave  from  8  would  enter  the  main  stream  right  on  the  tail  of 
the  wave  from  9,  the  arrival  times  at  1  would  be  ,855  and  1,  036. 

Thus,  all  arrival  times  in  this  pattern  would  be: 

Table  4 

Stamford,  Connecticut 
Summary  of  Automobile  Arrival  Time 

s  =  20  mph,  one  lane  traffic 


Origin 

Destination 

First  Car 

l.ast  Car 

7 

11 

.  782 

.832 

7 

12 

.946 

1.265 

10 

16 

.  170 

.  502 

10 

1 

.  146 

.  192 

9 

1 

.  524 

.855 

8 

1 

.855 

1. 036 
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B.  Foot  Traffic 


Calculation  of  v/alking  times  and  rates  is  a  much  simpler  proposition. 
Required  data  are  the  area  assignments,  and  the  maxinrim  and  minimum  dis¬ 
tances  for  each  origin-destination  combination.  The  distances  divided  by  the 
walking  rate  (here  taken  to  be  3  inph)  gives  the  time  of  arrival  of  the  first  and 
last  persons  from  each  origin.  Their  difference  is  the  arrival  interval  which, 
when  divided  by  the  number  of  persons  assigned,  gives  the  arrival  rate.  Table 
5  summarizes  these  data  for  Stamford,  Connecticut.  Independence  between 
automobile  and  pedestrian  traffic  is  assumed. 


Table  5 

Stamford,  Connecticut 
Inter-tract  Walking  Arrivals 


Dest, 

Origin 

Number 

Distance 
Min.  Max. 

Time* 

Start  Interval 

(persons /hr. ) 
Rate 

1 

13 

2987 

0 

1.  13 

0 

38 

7926. 90 

1 

15 

940 

0 

.  57 

0 

19 

4989.  24 

1 

21 

1600 

0 

1.  17 

0 

39 

4088. 80 

12 

8 

1355 

0 

1.  13 

0 

38 

3595.  90 

17 

21 

291 

0 

.  52 

0 

17 

1673.  25 

18 

20 

673 

0 

.  78 

0 

26 

2579.  81 

22 

21 

1933 

C 

1.  19 

0 

36 

5335.  08 

23 

16 

3411 

0 

.  82 

0 

27 

12387.  39 

3  mph,  in  1/100  hr.  units. 


C.  Combining  the  Data 

Perhaps  the  easiest  method  for  bringing  all  of  the  arrival  times  and 
rates  data  together  for  a  given  tract  is  to  draw  a  "bar  chart"  using  a  separate 
solid  line  to  indicate  on  a  time  scale  the  interval  duri.ig  which  traffic  is  arriv¬ 
ing  from  each  origin.  The  beginning  and  end  of  each  such  line  is  a  point  in 
time  at  which  the  arrival  rate  will  almost  certainly  change.  These  points  can 
very  easily  be  picked  off  such  a  chart,  and  the  arrival  rates  for  all  of  the  in¬ 
coming  waves  summed  to  give  a  total  arrival  rate  for  each  time  interval.  The 
end  result  would  then  appear  as  in  Table  b,  which  summarizes  the  results  of 
all  of  the  t  alculations  for  Census  Tract  1  in  Stamford,  Connecticut. 
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Table  6 


Stamford,  Connecticut 
Summary,  of  All  Arrivals  in  Tract  1 

20  mph  automobile  speed,  1  lane,  3  mph  walking  speed 


Destina¬ 

tion 

Period 

(hundredths  of  hrs.  ] 

1  Origins 

Arrivals 

(persons) 

Cumulative 

Arrival 

(persons) 

Arrival  Rate 
(persons/ 
hr. ) 

1 

0  - 

10 

13,  15,21 

1691 

1691 

17, 005 

1 

10  - 

15 

13,  15,21,  19 

1174 

2865 

23,  522 

1 

15  - 

19 

13,  15,21.  19,  10,5 

1412 

4277 

35, 798 

1 

19  - 

20 

13,21,  19.  10.5 

304 

4581 

30, 309 

1 

20  - 

38 

13,21,  19,5 

4509 

9090 

24.  952 

1 

38  - 

39 

21.  19,5 

172 

9262 

17, 025 

1 

39  - 

53 

19,  5 

1831 

11093 

12,936 

1 

53  - 

59 

19.  5,  9 

1141 

12234 

18,843 

1 

59  - 

79 

19,  9,  2 

3715 

15949 

18, 430 

1 

79  - 

82 

9.  2,  24 

560 

16509 

18,  546 

1 

82  - 

86 

9,  24,  6 

746 

17255 

18, 629 

1 

86  - 

104 

24,  8,  6 

3496 

20751 

19,418 

1 

104  - 

116 

24,  6 

1528 

22279 

12, 722 

1 

116  - 

144 

6 

1702 

23981 

6,  089 

1 

144  - 

184 

4 

2464 

26445 

6,  154 

1 

184  - 

226 

3 

2646 

29091 

6.  319 
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APPENDIX  4 


Glossaiy  of  Symbols 


Pedestrian  Traffic  Model 


area 


a  parametric  constant 

A^s  area  of  entire  community 

Ag  area  of  tracts  having  excess  shelter 

k  parametric  constant 

0  variable  of  integration 

P  population 

Pg  excess  population 

r  distance  from  center  of  core  shelter  area 

r^^  radius  of  core  shelter  area 

T2  radius  of  entire  community 

t  time 

V  walking  velocity 

Vehicular  Traffic  Model 

b  parametric  constant 

*^rrdn  least  distance  from  an  origin  to  a  destination  area 

H  number  of  households  in  origin  of  greatest  excess 

population  having  one  or  more  cars  av^ailable 

n  number  of  roads  leading  into  shelter  area 

Pg  excess  population 

Pill  excess  population  of  origin  with  greatest  excess  population 

Pmt  total  population  of  origin  with  greatest  excess  population 
Q  number  of  persons  per  vehicle 

R  unit  vehicle  arrival  rate 

S  vehicular  road  speed 

t  1 1  me 

t,  time  of  arrival  of  last  vehicle  on  artery  i 

t„-j  time  of  arrival  of  last  v'chic  le 

t,  time  of  arrival  of  first  vehic  le 
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APPENDIX  4 


1.  SUMMARY 

Two  independent  analytic  models  of  inter -tract  shelter  seeking 
traffic  are  presented  and  discussed:  one  for  pedestrian  and  one  for 
vehicular  traffic.  Methods  for  evaluating  the  model  parameters  are 
discussed.  Stamford  and  Waterbury,  Connecticut,  are  used  as  examples 
in  the  derivation  and  application  of  the  models.  The  results  are  com¬ 
pared  with  the  results  of  the  earlier  point -by -point  simulation  analysis  of 
these  two  cities. 


2.  INTRODUCTION  AND  DISCUSSION 
2.  1  Introduction 


The  analysis  of  the  flow  of  pedestrian  and  vehicular  traffic  enroute 
to  shelter  areas  was  a  basic  part  of  the  over-all  analysis  of  shelter-seeking 
operations  discussed  in  Chapter  II  of  this  report.  The  method  used  is  de¬ 
scribed  in  detail  in  Appendix  3.  It  is  essentially  a  paper -and -pencil  simu¬ 
lation  technique,  which  uses  actual  locations,  distances,  and  times  for  each 
group  of  people  and  shelters;  traces  the  flow  of  traffic  along  actual  routes 
on  a  point-by-point  basis;  and  resolves  conflicts  among  groups  by  means  of 
a  clearly  defined  set  of  rules  which  are  assumed  to  reflect  realistic  emer¬ 
gency  operating  conditions. 

While  simple  in  concept  and  techniques  of  execution,  simulation  is 
a  lengthy  and  tedious  procedure  better  suited  for  automatic  data  processing 
than  for  hand  computation.  Also  in  the  simple  form  used  here,  it  does  little 
to  reduce  the  problem  to  its  basic  elements  and  provide  any  real  understanding 
of  the  relationships  between  these  and  what  is  happening  in  the  real  world. 

For  these  reasons,  work  on  the  two  models  presented  in  this  appendix  was 
begun.  In  their  current  state  of  development,  these  models  already  point 
to  some  very  basic  relationships  between  community  parameters  and  the 
iiow  of  traffic  to  shelter  areas.  If  fully  developed  and  validated,  they  could 
provide  additional  insights  into  the  fundamental  processes  involved  in 
.she  Itc  r -seeking  operations,  and  could  be  used  to  estimate,  with  a  minimum 
of  (  alculation  and  data  collection,  the  flow  of  traffic  for  any  community. 


These  two  models  are  introduced  in  this  section  of  this  appendix  and 
the  results  obtained  using  them  are  compared  with  the  results  obtained 
using  the  previous  method.  The  pedestrian  traffic  model  is  then  for¬ 
mally  derived  in  Section  3  and  the  vehicular  traffic  model  in  Section  4. 
Section  5  contains  sample  calculations  for  Stamford  and  Waterbury, 
Connecticut. 

2.  2  Results 


Figures  1  through  8  on  the  following  pages  show  the  results  of  the 
work  discussed  in  this  appendix.  Figures  1  through  4  show  the  cumula¬ 
tive  arrival  curves  as  determined  by  both  the  analytic  method  developed 
herein  and  the  point -by -point  simulation  technique  previously  developed. 
Figures  5  through  8  show  population  density  and  rate  of  arrival  curves 
and  the  data  upon  which  they  are  based,  so  that  the  validity  of  the  basic 
assumptions  regarding  population  density  and  distribution  may  be  assessed. 

2.  2.  1  Figures  1  and  2:  Cumulative  Walking  Arii-vuls 

The  analytic  expression  developed  to  relate  walking  arrivals  to 

time  is: 


W(t)  =  P 

e 


r(ri  t 


for  all  values  of  a  not  equal  to  2.  In  the  special  cases  where  a  equals  2, 
the  function  takes  the  form: 


W(t) 


ln{  1  +  Vt/ r  j) 
'-nir^/rj) 


The  parameters  of  this  expression  are:  the  number  of  persons  who  must 
seek  shelter  outside  of  their  own  area  (Pg);  the  speed  at  which  these  people 
walk  (V);  a  measure  of  the  area  of  those  tracts  having  exce'-s  shelter  (rj); 
a  measure  of  the  area  of  the  entire  community  (r^);  and  a  parametric  con¬ 
stant  (W).  ^ 


^  Further  work  is  required  to  determine  whether  this  parameter  will  aiSume 
different  values  for  different  types  of  communities.  The  same  value;  <a~Z) 
was  used  herein  for  Stamford  and  Waterbury  with  equally  good  results 
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Figure  1.  Cumulative  Pedestrian  Arrivals  in  Shelter  Area 
vs.  Time- -Stamford,  Connecticut 


(sTiosjScj  JO  spuEsnoijx)  3/.tjBjnuin3 


Figure  2..  Cumulative  Pedestrian  Arrivals  in  Shelter  Area 
vs.  Time- -  Waterbury,  Connecticut 
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Stamford,  Connecticut,  Cumulative  Vehicular  Arrivals  in 
Shelter  Areas  vs.  Time:  20MPH  One  Lane  Traffic. 


.1 

1  1  1 

0,  5  1.0  1.4 

TIME  (HOURS) 

Figure  4. 

y/aterbur^,  Connecticut,  Cumulative  Vehicular 

Arrivals  in  Shelter  Areas  vs.  Time: 

20MPH  One  Lane:  Traffic. 
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ThJs  expression ’s  plotted  in  Figures  i  and  2  for  three  values  of 
walking  speed  for  each  community.  Also  plotted  is  the  cumulative  walking 
arrival  curve  determined  by  the  methods  of  Appendix  for  a  walking  speed 
of  3  mph.  (This  latter  curve  is  incomplete  in  that  it  does  not  extend  as  far 
in  time  as  the  first  three.  This  is  because  it  applies  only  to  those  persons 
who  have  one  mile  or  less  to  travel  to  their  snelter  areas.  Complete  com¬ 
parisons  between  the  two  methods  are  therefore  not  possible.  ) 

The  agreement  between  the  results  obtained  with  the  two  methods 
is  limited,  but  without  further  investigation  it  is  impossible  to  assert  with 
complete  confidence  that  either  is  more  valid  than  the  other.  We  do  be¬ 
lieve,  however,  that  the  discrepancies  between  the  two  methods  are  prin¬ 
cipally  due  to  the  inadequacies  of  the  earlier  method  and  that, despite 
obvious  limitat'ons,  the  method  desc  ribed  herein  will  prove  to  be  far  more 
useful  as  it  is  further  developed  and  applied. 

One  limitation  of  the  model  presented  here  is  that,  in  reality, 
people  will  still  be  arriving  from  the  fringes  of  the  community  after  the 
function  shows  100%  arrivals.  This  is  due  to  the  representation  of  the 
community  as  a  circular  area  (see  Section  3.  1).  Further  work  is  required 
to  determine  the  prec’se  effects  of  this  idealization  and  perhaps  to  investi¬ 
gate  the  utility  of  other  geometric  figures  as  approximations  to  the  true 
configuration  of  the  community.  Also,  the  results  should  be  validated 
against  the  actual  movement  of  people  rather  than  by  comparison  with  the 
first  crude  model  as  was  done  here. 

Addit’onal  work  which  would  greatly  increase  the  utility  of  the  model 
consists  in  the  determination  of  the  relationship  between  the  parameter  a  and 
known  characteristic  s  of  the  community  .  As  long  as  it  is  necessary  to  fit  the 
rate  cf  arrival  function  to  empir-cal  data,  the  use  of  this  model  will  require 
considerable  data  processing.  But  when  an  appropriate  value  of  a  can  be 
assumed  for  any  given  communHy,  it  will  be  possible  to  calculate  the  desired 
results  very  qu'ckJy  or  to  read  them  d’rectly  from  previously  prepared 
charts  and  graph®. 

The  curves  shown  in  Figures  1  and  2  show  that,  at  a  normal  walking 
speed  of  about  3  mpib,  it  is  possible  for  everyone  to  reach  shelter  in  about 
45  minutes.  They  also  show  that  the  number  of  arrivals  at  any  time  is  much 
more  influenced  by  delays  or  lowered  walking  speeds  than  by  any  feasible 
increase  :n  walking  speed.  We  can  therefore  tentatively  conclude  that  resi¬ 
dents  of  communities  such  a  s  Stamford  and  Waterbury  should: 

know  wh>^re  to  go.  The  distances  involved  are  small 
enough  that  everyone  can  reach  shelter  before  fallout 
a  r  r 'ves . 


4-8 


2.  get  moving  as  quickly  as  possible.  Time  lost  due  to 
delay  cannot  be  recovered  through  any  feasible  increase 
in  walking  speed. 

3.  walk,  don't  run.  The  value  of  additional  speed  is  slight, 
but  the  hazards  great.  Were  the  population  to  break 
into  a  run,  it  is  likely  that  fewer  people  would  safely 
reach  shelter  because  of  mob  casualties  and  disrupted 
traffic  flow. 

2.  2.  2  Figures  3  and  4:  Cumulative  Driving  Arrivals 

The  analvtic  expression  developed  to  relate  vehicular  arrivals  to 

time  is: 


D(t)  = 


P 

e 

1  +  nb 


where  is  the  number  of  persons  who  must  use  vehicles  to  travel  to 
shelter  areas;  n  is  the  number  of  main  traffic  arteries  leading  into  shelter 
areas;  b  is  a  measure  of  population  groupings  in  the  community;  t^  is  the 
time  of  arrival  of  the  first  vehicle;  and  tm  the  time  of  arrival  of  the  last 
vehicle. 


Perhaps  the  outstanding  feature  of  this  model  is  the  existence  of  an 
explicit  relationship  between  the  parameter  b  and  observable,  measurable 
characteristics  of  the  community  (see  Section  4.  3.  1).  Thus,  it  is  possible 
to  apply  this  model  directly  to  any  community  without  having  to  fit  the  rate 
of  arrival  curve  to  desc^^iptive  community  data  as  is  currently  necessary 
with  the  pedestrian  traffic  model. 

The  vehicular  traffic  flow  model  is  applied  in  Figures  3  and  4  to 
Stamford  and  Waterbury  for  those  areas  determined  to  be  origins  and 
destinations  in  Section  of  the  report.  The  areas  and  the  total  numbers 
of  people  and  shelters  involved  are  therefore  not  identical  to  those  used 
in  applying  the  pedestrian  traffic  flow  model.  The  difference  is  the  exclu¬ 
sion  here  of  all  persons  residing  within  origin  tracts  less  than  one  mile 
center-to-center  distant  from  their  destination.  This  was  done  so  that 
the  results  here  would  be  comparable  with  the  results  of  the  earlier  vehicu¬ 
lar  traffic  flow  analysis  results,  which  are  as  close  an  approximation  of 
reality  as  is  available  to  us  at  this  time. 
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The  cumulative  arrival  curves  are  plotted  for  a  driving  speed  of 
20  mph.  The  form  of  both  curves  compares  very  favorably  with  that  from 
our  previous,  more  detailed  analysis.  The  difference  between  the  current 
and  previously  calculated  function  for  Waterbury  is  due  primarily  to  the 
differences  in  the  estimates  of  the  time  of  arrival  of  the  last  vehicle. 

Were  there  greater  precision  of  significant  value  in  a  particular 
case,  it  might  be  worthwhile  to  employ  the  method  of  Appendix  to  deter¬ 
mine  the  cumulative  vehicular  arrival  function.  The  difference  in  the  time 
required;  however,  in  general  will  give  the  current  model  far  greater 
utility.  Once  the  assignments  have  been  made,  the  main  arteries  selected, 
and  the  v  alues  of  the  constants  determined,  the  entire  function  can  be  calcu¬ 
lated  and  plotted  by  hand  in  an  hour  or  twoJand  previously  prepared  para¬ 
metric  curves  could  greatly  reduce  calculation  time.  The  method 
of  Appendix  involves  over  a  week  of  hand  computation  for  a  city  the  size 
of  these  two. 

No  new  conclusions  have  evolved  from  this  effort  relative  to  vehicu¬ 
lar  traffic  flow,  both  because  this  model  agrees  so  well  with  the  previous 
model,  and  because  there  has  been  no  opportunity  to  take  advantage  of  its 
analytic  structure  to  determine  the  effects  of  systematically  varying  values 
of  the  community  parameters  related  to  vehicular  traffic  flow. 

2.  2.  3  Figures  5  and  6:  Population  Density  Curves 

The  basic  postulate  of  the  pedestrian  traffic  model  is  that  the  density 
of  the  excess  population— those  persons  who  must  seek  shelter  outside  of 
their  own  tract  area  — is  inversely  related  to  the  distance  from  the  shelter 
area  to  which  they  would  go.  To  the  extent  that  this  postulate  is  valid,  our 
confidence  in  the  validity  of  the  model  is  increased;  and  to  the  extent  that 
this  postulate  is  not  descriptive  of  reality^  our  confidence  in  the  validity  of 
the  model  js  lessened. 

In  Figures  5  and  6-  the  postulate,d  relationship  is  given  by  a  smooth 
curve  fitted  to  the  data  by  using  the  values  of  the  parameters  computed  in 
Section  5.  1  and  5.  2.  The  actual  community  data  are  given  by  the  plotted 
points.  Straight  lines  connecting  the  points  have  been  drawn  as  an  aid  to 
visualiza,,  on.  Eac.h  pomt  plots  the  distance  from  the  center  of  an  origin 
tract  to  the  ceiiicr  )f  the  nearest  shelter  tract  against  the  excess  population 
of  that  origin  tract.  Where  no  one  shelter  U'-ea  can  accept  all  of  the  excess 
population  of  a  given  origin,  multiple  points  appear.  While  these  curves 
might  or  might  not  represent  a  "best"  fit.  we  feel  that  the  degree  of  fit  is 
.startling,  and  does  give  one  confidence  in  the  over-all  validity  of  the  model. 
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Distance  From  Center  of  Core  Shelter  Area  ^Miles) 

Figure  5  Population  Density  vs.  Distance  From  Center  of 
Core  Shelter  Area,  Stamford,  Connecticut  . 
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Figure  6.  Population  Density  vs.  Distance  from  Center  of  Core 
Shelter  Aiea.  Waterbury,  Connecticut. 


2.  2.  4  Figures  7  and  8:  Rate  of  Arrival  Curves 


The  basic  postulate  in  the  vehicular  traffic  model  is  analogous  to  the 
population  density  postulate:  we  assume  regularity  in  the  difference  in  the 
number  of  vehicles  originating  at  each  of  the  centers  of  excess  population. 

That  regularity  is  described  by  the  postulated  mathematical  function. 

Figures  7  and  8  show  the  postulated  functions  for  Stamford  and 
Waterbury  as  continuous,  smooth  curves,  and  the  actual  number  of  vehicles 
arriving  at  the  shelter  areas  as  a  step  function.  The  step  function  was 
generated  by  plotting  the  time  period  over  which  automobiles  would  be  arriv¬ 
ing  over  each  of  the  main  arteries,  assuming  simultaneous  initial  arrivals. 

Note  that,  while  the  curves  do  indeed  fit  well,  the  meaning  is  not 
the  same  here  as  w'ith  the  pedestrian  traffic  model.  There,  the  data  fitted 
were  basic  community  data,  and  the  implication  was  that  the  model  was  a 
valid  one.  Here  the  data  fitted  are  derived  using  the  methods  of  Appendix 
3  ,  and  the  implication  is  that  the  model  developed  herein  is  a  valid  approx¬ 
imation  to  those  methods.  Any  inference  about  the  basic  validity  of  this 
model  must  be  based  either  on  an  acceptance  of  the  method  of  Appendix  3 
or  upon  empirical  validation  procedures. 

3  FORMAL  DEVELOPMENT  OF  THE  PEDESTRIAN  TRAFFIC  MODEL 
3.  1  The  Idealized  Community 

For  computational  convenience,  the  community  is  represented  by  a 
pair  of  concentric  circles  with  radii  rj^irj  .  r^^  is  selected  so  that  the 

area  of  the  circle  it  defines  is  equal  to  the  area  of  thosv^  tracts  with  an  excess 
of  shelters.  r^  is  selected  so  that  it  defines  an  area  equal  to  the  area  of  the 
entire  community.  Figures  9  and  10  show  the  relationship  between  these 
idealized  areas  and  the  actual  areas  for  Stamford  and  Waterbury,  Connecticut. 

Intuitively,  it  would  appear  that  t’lis  geometric  approximation  is  a 
severe  departure  from  reality,  particularly  in  the  case  of  Stamford.  Exam¬ 
ination  of  Figures  5  and  6  in  the  previous  section,  however,  shows  that  it 
leads  in  both  cases  to  a  rather  good  approximation  of  tiiat  characteristic 
of  the  community  which  we  are  most  interested  in  approximating  — the  density 
of  the  population  at  any  given  distance  from  the  center  of  the  core  shelter 
a  ’•ea . 
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P'igure  10.  Watcrbury,  Connecticut.  Actual  and  ideali*ed 
. I  road  of  excesti  population  and  shelter. 
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3.  2  Mathematical  Development 


If  we  assume  that  the  density  of  the  excess  population  in  the  area 
between  and  rg  decreases  with  increasing  distance  from  the  center  of 
the  core  shelter  area  in  accordance  with 


where  k  and  a  are  parametric  constants,  then  we  can  make  the  following 
substitution  to  express  equation  (1)  in  polar  coordinates 


dA  =  rdrd0 


and  solve  equation  ( 1)  for  dP 


dP  =  kr^'^drdg 


from  which  it  follows  that 


=  k 


4rd0 


which,  when  integrated,  becomes 


P  =  2,rk  +  C 

2-a 


a  #  2 


Equation  (5)  describes  a  family  of  curves  parametric  in  a.  Each 
of  these  curves  must  meet  two  real-world  conditions  which  permit  the 
evaluation  of  k  and  c.  The  first  of  these  conditions  is  that  none  of  the 
excess  population  is  contained  within  the  inner  circle  which  defines  the 
"core  shelter  area"  of  the  community.  That  is, 


Condition  1: 


Plr^)  =  0 


'rom  equmon  (S); 


P(r,)  - 


2-a 


t  C 


Sol  VI !  g  for  C: 


-  2ffk  r. 
2-a 
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a  ¥=  2 


(6) 


and  substituting  in  equation  (5): 


P  = 


2ffk 

2-a 


(r 


a-a. 


rj-a 


) 


Now  k  is  evaluated  by  considering  the  second  condition— that  all  of 
the  excess  population  is  contained  within  the  outer  circle  which  defines 
the  limits  of  the  community. 

Condition  2:  =  P^ 


from  equation  (6): 


2ffk  .  j.a  a“a\ 

P(ra)  =  (rf  -  rf  *) 


Solving  for  k: 


k  = 


Pe(2-a) 


and  substituting  in  equation  (6): 


(7) 


P  = 


r?~*) 

(r^^-rj"^) 


a  #  2 


(8) 


Equation  (7)  gives  the  excess  population  contained  within  some  circle 
defined  by  radius  r,  where  r^  ^  r  rj  .  We  arc  interested  in  the  time  re¬ 
quired  for  this  population  P  to  reach  the  boundary  of  the  core  shelter  area. 
Therefore,  since 


t= 


(9) 


where  V  is  the  walking  velocity,  we  substitute 


r  ^  Ti  f  Vt 

in  equation  (7),  which  yields  the  desired  result: 


W(t) 


a  #  2 


(10) 
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For  the  special  case  a=  2,  the  following  result  may  be  derived  by 
taking  the  limit  of  equation  (10)  as  a  approaches  2; 


W(t)  =  Pg 


ln(i  Vt/r,) 
Intrg'rj) 


a  =  2 


(11) 


This  form  of  tlie  expression,  equation  (11),  was  found  to  give  a  good 
fit  to  the  data  for  both  Stamford  and  Waterbury,  and  was  therefore  used  in 
all  calculations  in  Section  5  of  this  appendix. 

3.  3  Evaluating  the  Parameters 

3.  3.  1  a:  The  Parametric  Constant 


The  appropriate  value  of  a  must  be  determined  by  standard  curve¬ 
fitting  procedures,  as  no  relationship  between  it  and  other  parameters,  such 
as  that  found  for  b  in  Section  4.  3.  1,  has  yet  been  found.  Should  the  value 
a  =  2  yield  a  good  fit,  as  is  the  case  for  the  two  cities  discussed  in  this  report, 
then  the  computationally  simpler  form  of  the  final  expression  given  as  equation 
(11)  may  be  used. 

3.  3.  2  ri:  Radius  of  the  Core  Shelter  Area 

The  radius  of  the  core  shelter  area  is  determined  by; 

ri  =  /  A  '/ff-  ( 1 2) 

9 

where  is  the  sum  of  those  areas  having  excess  shelter  spaces. 

3.  3.  3  fa:  Radius  of  the  Entire  Community 

The  radius  of  the  entire  community  is  given  by: 

r,  ^  /imr  (13) 

"  c 

where  A^.  is  the  total  area  of  the  community,  including  areas  both  of  excess 
shelter  and  of  excess  population. 

3.  3.  4  V:  Walking  Velocity 

The  value  of  walking  velocity  is  chosen  as  a  planning  parameter.  In 
this  appendix,  walking  arrival  curves  are  given  for  one,  three,  and  five 
miles  per  hour  walking  speed,  and  compared  with  the  results  of  our  earlier 
analysis  which  assumed  three  miles  per  hour. 
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4.  FORMAL  DEVELOPMENT  OF  THE  VEHICULAR 

TRAFFIC  MODEL 

4.  1  The  Idealized  Community 

In  any  community,  the  bulk  of  the  vehicular  traffic  is  effectively 
restricted  under  normal  conditions  to  a  relatively  few  main  arteries.  We 
assume  that  this  will  be  the  case  also  during  a  civil  defense  emergency. 
Figure  11  shows  the  flow  of  vehicular  shelter- seeking  traffic  for  Stamford, 
Connecticut,  based  upon  this  same  assumption  in  our  earlier  analysis. 

Note  that,  although  many  roads  are  used,  the  traffic  rapidly  converges  on 
only  eight  arteries  which  actually  enter  the  shelter  areas. 

Fig  ure  12  shows  the  idealized  community  configuration  upon  which 
the  vehicular  traffic  flow  model  is  based.  Each  group  of  tracts  having  insuf¬ 
ficient  shelter  and  using  the  same  main  artery  to  enter  a  shelter  area  is  con¬ 
ceived  of  as  a  single  origin  reaching  the  shelter  area  by  a  single  connecting 
route.  The  population  of  each  origin  is  the  sum  of  the  populations  of  these 
tracts.  All  of  the  tracts  with  an  excess  of  shelter  areas  are  conceived  of 
as  a  single  "core  shelter  area.  "  The  differences  in  the  length  of  the  dif¬ 
ferent  connecting  routes  are  ignored,  as  they  have  proved  in  practice  to  be 
negligibly  small.  Therefore,  vehicles  from  all  of  the  origins  will  begin  to 
arrive  simultaneously  in  the  core  shelter  area. 

4.  2  Mathematical  Development 

Fig  ure  13  shov's  the  discrete  function  describing  the  assumed  traffic 
flow.  If  the  vehicles  which  begin  to  leave  each  of  the  origins  simultaneously 
at  time  zero  and  at  rate  R,  the  rate  of  arrival  at  the  destination  is  zero  until 
time  tjj,  at  which  time  vehicles  begin  to  arrive  simultaneously  on  all  main 
arteries.  If  there  are  n  such  arteries,  Q  persons  per  car,  and  the  vehicle 
flow  rate  is  R  cars  per  artery  per  unit  »‘rne,  then  the  total  initial  rate  of 
arrival  is  nQR  persons  per  hour.  This  .ate  of  arrival  holds  until  all  of  the 
people  have  arrived  from  the  origin  of  least  excess  population,  at  which  time 
the  rate  of  arrival  drops  discretely  to(n-  1)QR.  This  process  continues 
until  time  t,^^  where  the  last  car  from  the  origin  with  the  greatest  excess 
population  arrives,  and  the  rate  of  arrival  drops  dis*  rrtely  from  UR  to  zero. 

Th«  (o*a1  nuii'ber  of  persons  involved  *s  given  hy: 

~  UR  £tt.  -  i  )  t  1.  i . n  (14) 

I  o 

where  1.  is  the  duration  of  traffic  flow  on  the  i-lh  main  artery.  In  the  special 

1 

case  where  the  traffic  flow  on  succeeding  arteries  is  exhausted  at  regular 
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Figure  13.  Discrete  Rate  of  Arrival  Function 


intervals,  traffic  continues  to  flow  on  each  artery  fcr  exactly  (l/n)(tj^-  t^) 
time  units  longer  than  on  the  preceding  one  before  becoming  exhausted. 
Therefore 

t.-t  =  -i2— 2-  (n-j)  j=n+!-i  (IS) 

1  o  n 


and 


n 


n 

r  (n-j) 

J=1 


QR{t  - 
m 


(H-n) 

2 


(16) 


The  discrete  rate  of  arrival  function  discussed  above  can  be  approxi¬ 
mated  by  the  continuous  function  (see  Figure  14): 

t  - 1  .b 


D'(t)  =  OR  [(1  -  n)  t  n] 

in 


t  <  t  :S  t 
o  m 


(17) 


where 

t  =  time 

Q  =  number  of  persons  per  vehicle 

R  =  vehicle  arrival  rate  in  cars  per  artery  per  unit  time 

n  =  number  of  arteries  leading  into  shelter  areas 

t  =  time  o*  arrival  of  first  vehicle 
o 

t  =  time  of  arrival  of  last  vehicle 

m 

b  =  a  parametric  constant 

which  has  the  following  properties 


D'(t  )  =  nOR  (18.  1) 

o 

D'(t  )  =  OR  (18.2) 

m 

D(t  )  =  P.  (18.3) 

m  ® 

These  properties  may  be  verified  by  inspection.  The  quantity 
(t-  to)/(tm~  to)  will  vax-y  from  zero  to  one  as  t  varies  from  to  to  tm* 

When  it  is  equal  to  zero,  the  value  of  the  entire  expression  is  nQR;  when 
it  is  equal  to  one,  the  value  of  the  entire  expression  is  OR.  The  third 
property  (18.  3)  is  used  to  fit  the  curve  to  the  data  for  any  given  community. 
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Examination  of  the  geometry  of  Figure  13  will  make  clear  the  relationship 
between  D(t)  and  Pg.  In  the  special  case  where  the  traffic  flow  on  succeeding 
arteries  is  exhausted  at  regular  intervals,  the  parameter  b  is  equal  to  unity 
and  the  area  under  the  linear  approximation  is  equal  to  tg)  (l+n)/2 

which  is  in  agreement  with  equation  (16). 

In  general,  however,  the  function  is  not  linear  and  it  is  necessary  to  use 
empirical  data  to  calculate  the  value  of  b  (see  Sections  5.  3  ai-d  5.  4  of  this  appendix. 
When  b  is  less  than  unity,  the  function  is  concave  upward;  when  b  is  greater 
than  unity,  the  function  is  concave  downward.  These  values  of  the  parameter 
correspond  respectively  to  the  real-world  cases  where  the  total  number  of 
vehicles  is  less  evenly  distributed  and  more  evenly  distributed  among  the 
several  main  arteries. 


The  cumulative  arrival  function,  which  is  what  we  are  really  interested 
in,  is  obtained  by  integrating  D'(t)  dt  subject  to  the  initial  condition  expres¬ 
sed  in  (18.  3). 

p  t-t  b 

D(t)  =  OR  [U  -  n)  (^— ^  }  +  nj  dt  .  (19) 


First  the  appropriate  substitutions  are  made  in  order  to  integrate  with 


Jjirst  the  appropriate  substi 
respect  to  I  •'  t-t 


t  - 1 
m  o 


t-t  ,b 


D(t)=QR(t^-to)J  [(l-n){^}  ^n]-^ 

w  mo  Tin 


m  o 


and  the  expression  is  integrated: 


D(t)  =  - 


QR(t  -t  ) 
m  o 


t-t  b+ 1 


m  o 


m  o 


then  property  ( 183)  is  employed: 


P  =  D(t  )  =  QR(t  -t  )-71T^ 
e  m  m  o  ( 1+  b; 
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4.  3.  3  R:  Vehicle  Arrival  Rate 


The  vehicle  arrival  rate  for  each  artery  is  based  upon  an  estimate 
of  the  capacity  of  the  roads  leading  into  the  shelter  areas.  The  arrival  rate 
for  single  lane,  20  niph  traffic  (1650  vph)  derived  earlier  is  used  here. 

4.3.4  n:  Number  Arteries  Leading  Into  Shelter  Area 

The  number  of  arteries  leading  into  the  shelter  area  is  determined 
by  studying  the  area,  and  deciding  which  roads  as  a  matter  of  policy  will  be 
designated  main  arteries.  In  many  cases,  local  geography  and  street  lay¬ 
out  will  limit  the  choices  available;  in  all  cases,  it  will  be  necessary  to 
consider  traffic  and  other  planning  activities. 

4.  3.  5  tot  Time  of  Arrival  of  First  Vehicle 


to  is  determined  by  the  relationship: 


mm 


(28) 


where  d^^^  is  the  least  distance  from  an  origin  to  a  destination  area  and  S 
is  the  road  speed  of  the  vehicles. 

4.  3.  6  tm*  Time  of  Arrival  of  LASt  Vehicle 

tin  determined  by  the  relationship: 


_  X  ni 

^m  m  QR 


(29) 


where  dm  c..  stance  from  the  origin  of  greatest  excess  population  and 

Pm  is  the  excess  population  of  that  origin. 

5.  COMPUTATIONS 


5.  1  Stamford  Pedestrian  Traffic 

The  total  shelter  area  in  Stamford,  Connecticut,  comprising  eight 
of  the  24  census  tracts  in  that  city,  is  4.52  square  miles.  The  area  of  the 
entire  community  is  37.  53  square  miles.  ^  Therefore  using  equations  (12)  and  (13): 


^Areas  were  obtained  by  direct  measurement  of  the  maps  appearing  the  PHC 
series  publications  of  the  1960  Census  of  the  U.  S.  and  may  differ  slightly  from 
those  published  elsewhere.  This  does  not  affect  the  validity  of  the  results  obtained. 
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=  i/4.  ft27fr  =  1.235  miles 
Tg  =i/iy.  n'  -  3.  455  miles 


Normally,  it  is  not  necessary  to  compute  the  value  of  k.  It  is 
done  here  in  order  that  the  relationship  assumed  in  equation  (1)  can  be 
plotted  to  demonstrate  the  degree  of  fit  realized  by  assuming  the  value 
of  a  equal  to  2.  The  limit  of  equation  (7)  as  a  approaches  2  (the  value 
assumed  herein),  is  given  by: 

P 

k=:  ^ 

2jrln{r^/rg) 

The  excess  population  in  Stamford  numbers  44,  397.  Therefore 


44397 
27r(l.  03) 


6,850 


The  expression  for  the  rate  of  arrival  is  therefore 

dP  =s  685(*r"2 

and,  using  equation  (11),  the  expression  for  the  cumulative  number  of  arrivals  is 

P=  430001n(l+~^) 

The  rate  of  arrival  curve  appears  as  Figure  5  on  page  4-11.  The  cumula¬ 
tive  arrival  curves  are  plotted  for  three  values  of  V  in  Figure  1  on  page 
4-4. 


5.  2  Waterbury  Pedestrian  Traffic 

The  shelter  area  in  Waterbury  is  5.  46  square  miles,  and  the  area 
of  the  entire  community  is  30.  12  square  miles. ^  Therefore 

~  /5.  46/Tr'  =  1.  32  miles 
Tg  =/30.  12  /  TT  =  3.  10  miles 

^See  footnote  on  previous  page. 
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Again  assuming  a=  2,  using  equation  (30)  and  substituting  the 
Pg  =  68,  042, 


68042 

2ff{.851) 


12,  700 


dP  =  12,  700r“* 


and,  using  equation  (1 1), 

P=r  80,000Xn(l+-Y^) 


These  expressions  are  plotted  in  Figure  6,  page  4-12,  and  in  Figure  2, 
page  4-5. 


5.  3  Stamford  Vehicular  Traffic 

From  the  I960  Census  of  the  United  States,  the  National  Shelter 
Survey  Phase  II  results,  and  a  detailed  map  obtained  from  the  city 
engineering  department,  the  following  values  have  been  obtained;^ 


Pe 

31,  307 

P 

m 

13,  138 

P 

mt 

= 

13,  569 

d  . 
min 

.  130  miles 

d 

m 

.  306  miles 

n 

= 

8 

H 

SZ 

3.  586 

^Implicit  in  these  figures  is  an  assignment  of  persons  from  origins  to 
certain  destinations  and  to  movement  on  foot  or  by  vehicle.  The  assign¬ 
ment  made  here  is  the  same  as  in  our  earlier  analysis  so  that  comparisons 
can  be  made  between  the  two  approaches.  This  i.«  why  the  value  of  P^  here 
differs  from  th^  ..  used  in  the  application  of  the  walking  model. 
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From  these,  the  following  were  calculated: 


13,  138  »  31,  207 
31,  207  -  8(13,  138) 


.  2445 


t 

o 


S 


1.30 

20 


.  065  Hrs. 


3.06  13,  138 

20  (3.7839)(1650) 


2.  257  Hrs. 


13,  569 
3,  586 


3.  7839 


The  rate  of  arrival  function  for  Stamford  is  determined  by  substituting  the 
above  in  equation  (17): 


D^(t)  =  (3.  7839) (1650)  |(i-8){- 


t  - .  065 


2.  257-  .065 


}■ 


2445 


065 <  t<  2.  257 


D^(t)  =  6243 


t  -  .  065 
2.  192 


.  2445 


SJ 


.  065t  <  2.  257 


and  the  cumulative  arrival  function  is  similarly  determined  by  substitution 
in  equation  (23); 


Djlt) 


=  31,  207 


[(1  -  8)  {■ 


t  -  .  065 
2.  257 -.065 


1  +  .  2445 


+  8(1 +.  2445)-^^ 


t-.065 


257-.  065/ J 


l+8(.2445) 


Dg(t)  =  1056 


r  t  -  . 065 
i.  2.  192 


1.  2445 


t-.  065  I  " 

2.  192  J  J 


.  065t  <  2.  257 
.  065  S.t  ^2.  257 


These  curves— the  rate  of  arrival  and  the  cumulative  arrivals— appear  as 
Figure  7,  page  4>14,  and  Figure  3,  page  4-6. 
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5. 4  Waterbary  Vehicalar  Traffic 

From  the  same  sources  as  for  Stamford,  the  following  data  were 
obtained^ 


p 

e 

33,  401 

P 

m 

7,  698 

**mt 

■  = 

8,482 

*^min 

s 

1.  30  miles 

d 

m 

— 

1.  70  miles 

n 

10 

H 

Si 

2,  088 

From  these,  the  following  were  calculated: 


7698-  33401 
33401-  10(7698) 


.  5898 


_  ^min  _  1. 30 
^o”  S  ”20 


.065  Hrs. 


JSl  - 

■  H  ”  2088 


4.  062 


1.70  7698 

20  (4.062)(1650) 


1.  234  Hrs. 


^See  footnote  on  Page  4-32. 
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The  rate  of  arrival  curve  for  Waterbary  ie,  by  eubstitution,  in  equation 
(IT): 


oy  t)  =  (4.  062)11650)  [(1  -  10)  { 


.  5898 

}  ^  ‘“J 


.  065  <  t  <  1.  234 


D^(t)  =  6702  [9  { 


t-.065  Y 
1.168  / 


5898 


+  10  J 


.  065  <  t  <  1.  234 


and  the  cumulative  arrival  function  ia,  by  sub8titution>  equation  (23): 


D^(t)  ~ 


~  33401 


[a-i0)(_ 


t  -  .065 


1  +  .  5898 


}  .  10,1..  5898 


234-.  065 


1+  10  (.5898) 


.  065  <  t  <  1.  234 


4842 


t-065  **’« 

1.  169  J 


+ 


15.898  -[ 


t-.  065  1  1 
1.  169  J  J 


.  065  <  t  <  1.  234 


Theae  curvea  appear  aa  Figure  B,  page  4>15,  and  aa  Figure  4.  page  4>9. 


4.35 


APPENDIX  5 


APPENDIX  5 


Model  to  Predict  Niimber  of  Tract  Residents  Sheltered  Within 
Tract  as  a  Function  of  Elapsed  Time  From 
When  Residents  Leave  Homes 


This  model  will  be  developed  with  the  aid  of  an  illustration.  The  illus¬ 
tration  is  for  census  tract  #14  of  the  City  of  Stamford,  Connecticut.  The 
basic  data  required  by  the  model  is  as  follows: 

Residents  to  be  sheltered  in  the  tract:  4781 

Nximber  of  shelters  in  tract:  10 

Tract  area:  .  642  mi.^ 

2 

Population  density  of  tract:  7450  people  per  mi. 

It  is  assumed  that  the  number  of  residents  per  unit  area  is  constant 
throughout  the  tract,  and  it  is  also  assumed  that  the  shelters  are  so  spaced 
that  an  area  and  hence  a  fixed  number' of  residents  can  be  associated  with 
each  shelter.  For  the  present  illustration  the  sub-area  associated  with  each 
shelter  is  .  0642  sq.  miles,  and  the  number  of  residents  in  this  sub-area 
is  478.  1.  It  is  also  assumed  that  the  residents  within  a  shelter  sub-area 
will  always  proceed  to  that  shelter  first  and,  if  it  is  full,  will  then  proceed 
to  the  next  nearest  shelter. 

It  is  convenient  to  rank  the  shelters  in  order  of  their  capacity  as  shown 
in  Figure  1.  When  the  warning  is  received  by  the  residents,  they  will  pro¬ 
ceed  to  their  nearest  shelter.  Thus,  478  residents  will  "visit"  each  of  the 
ten  tract  shelters  and,  as  can  be  seen  from  Figure  1,  some  of  the  residents 
will  find  a  shelter  space  available  while  others  will  not.  Table  1  shows  that 
3,  474  residents  can  expect  to  find  shelter,  while  1,  307  will  have  to  visit  at 
least  one  more  shelter  to  find  space.  The  situation  faced  by  the  1,  307  un¬ 
sheltered  residents  is  that  three  shelters  are  not  full  and  seven  shelters  are 
full.  However,  each  resident  has  visited  a  full  shelter  and  presumably  is 
aware  of  its  location  so  will  proceed  to  one  of  the  remaining  nine  shelters. 
Because  six  of  the  lemaining  shelters  are  full,  6/9  of  the  unsheltered  resi- 
dentSj  i.  e. ,  871,  will  be  unsheltered, after  these  second  visits  to  a  shelter, 
Three-nonths  of  the  unsheltered  residents,  i.  e. ,  436,  will  visit  shelters 
vthich  are  not  completely  full.  It  ^  assumed  that  these  436  residents  will 
divide  evenly  among  the  three  partially  full  shelters,  or  145  residents  per 
shelter.  Returning  again  to  Figure  1,  it  can  be  seen  that  145  residents  per 
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Figure  1.  Rank  of  Shelters  in  Order  of  Capacity- - 
Tract  No.  14. 
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partially  filled  shelter  fills  up  one  mor«  shelter  (Rank  8)  and  leaves  87 
residents  more  still  unsheltered.  Thus,  there  are  871  plus  87  or  958 
residents  that  must  make  a  third  visit  to  a  shelter.  The  situation  faced 
by  these  958  residents  is  two  shelters  partially  filled  and  8  shelters  com¬ 
pletely  full.  However,  each  resident  has  now  visited  two  full  shelters  and 
presumably  will  not  return  to  them.  Thus,  each  resident  now  has  a  "popu¬ 
lation"  of  eight  shelters  to  visit,  two  of  which  are  partially  full  and  six  of 
which  are  completely  full.  Again,  it  is  assumed  that  6/8  (i.  e.  ,  718)  of  the 
residents  visit  full  shelters  and  must  make  a  fourth  visit.  240  residents 
visit  shelters  which  are  partially  empty,  or  120  residents  per  shelter.  It 
can  be  seen  from  Figure  1  that  120  residents  do  not  fill  up  either  of  the 
shelters,  hence  all  240  are  now  sheltered,  leaving  718  to  make  a  fourth 
visit.  - 

The  procedure  is  now  continued  until  all  residents  are  sheltered.  The 
number  of  residents  sheltered  after  each  visit  in  contained  in  Table  1,  which 
follows: 


Table  1 

Cumulative  Residents  Sheltered 


Shelter 

Visits 

Residents 
Unsheltered 
Before  Shelter 
Visit 

Residents 
Sheltered 
After  Visit 

Residents 
Requiring 
Another 
Shelter  Visit 

Shelters 
Partially 
Full  at  end 
of  Visit 

Shelters 
Completely 
Full  at  end 

of  Visit 

1 

4781 

3474  ^ 

1307 

3 

7 

2 

1307 

349 

958 

2 

.8 

3 

958 

240 

718 

2 

8 

4 

718 

129 

589 

1 

9 

5 

589 

98 

491 

1 

9 

6 

491 

98 

393 

1 

9 

7 

393 

98 

295 

1 

9 

8 

295 

98 

197 

1 

9 

9 

197 

98 

99 

1 

9 

10 

99 

99  . 

0 

0 

10 

The  next  step  in  the  model  is  to  compute  the  number  of  residents  sheltered 
as  a  function  of  time.  For  purposes  of  this  computation,  it  is  assumed  that  the 
sub-  area  associated  with  a  shelter  is  in  the  form  of  a  circle.  Each  sub-area 
contains  .  0642  mi.^,  which  corresponds  to  a  circle  of  .  143-mile  radius.  Return¬ 
ing  now  to  Figure  1,  it  can  be  seen  that  130  residents  visit  each  shelter  before 
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the  first  shelter  is  full.  From  the  population  density  of  the  tract,  it  is  found 
that  these  residents  come  from  an  area  of  .  0174  sq.  miles  surroimding  the 
shelter.  The  radius  of  this  area  is  .  0745  miles.  Thus,  the  first  of  this  group 
of  130  residents  travels  zero  miles  to  the  shelter  while  the  last  travels  .  0745 
miles.  At  a  speed  of  3  mph,  elapsed  time  for  the  last  of  the  group  to  be  shel¬ 
tered  is  .0247  hours. ^  Because  there  are  ten  shelters,  1300  residents  are 
sheltered  in  the  time  period  from  zero  to  .  0247  hours  after  warning.  A 
second  shelter  is  filled  after  87  more  residents  visit  each  shelter.  Again 
from  population  density,  it  is  foxind  that  this  group  of  residents  may  have  to 
walk  as  far  as  .  0962  miles  to  get  to  the  shelters.  There  are  now  nine  shel¬ 
ters  which  accept  these  residents  for  a  total  of  783  people.  At  3  mph,  the 
time  period  during  which  this  group  of  783  people  is  sheltered  extends  from 
.  0247  hours  to  .  0321  hours.  Thus,  at  the  end  of  .  032  hours  the  total  number 
of  residents  sheltered  is  1300  +  785,  or  2083.  This  process  is  continued 
tmtil  the  nximber  of  residents  finding  shelter  on  the  first  shelter  visit  and 
the  time  period  during  which  they  find  shelter  are  determined. 

It  is  now  necessary  to  determine  the  time  required  to  find  shelter  for 
those  residents  who  are  not  sheltered  on  their  first  visit  to  a  shelter.  There 
are  478  residents  associated  with  each  shelter  in  the  census  tract  and,  as 
was  noted  earlier,  one  shelter  is  filled  by  130  residents  leaving  348  residents 
to  visit  this  shelter  and  then  proceed  to  the  next  nearest  shelter.  The  time 
required  for  the  shelter  to  fill  is  .  0247  hours,  and  .  0230  hours  later  (i.  e. , 

.  0477  hours,  from  departing  homes)  each  resident  has  visited  the  shelter 
with  which  he  has  been  associated.  The  time  interval  of  .  0230  hours  also 
represents  the  interval  during  which  residents  who  have  not  found  space  in 
the  first  shelter  visited  are  arriving  at  a  second  shelter. 

The  minimxim  distance  between  shelters  in  this  illustration  is..  286  miles, 
therefore,  the  minimum  distance  traversed  in  visiting  the  first  two  shelters 
is  .  360  miles  and  the  maximum  distance  travelled  is  .429  miles.  At  3  mph, 
these  distances  correspond  to  .  120  hours  and  .  143  hours  respectively,  which 
means  that  all  residents  who  find  shelter  space  on  their  second  visit  to  the 
shelter  are  sheltered  in  the  time  interval  between  .  120  hours  and  .  143  hours 
after  departing  their  homes.  From  Table  1  it  can  be  seen  that  349  residents 
are  sheltered  during  this  time  period.  The  time  period  during  which  residents 
find  shelter  space  during  their  third  visit  to  a  shelter  is  determined  by  adding 
the  transit  time  between  shelters  (i,  e. ,  .  095  hours)  to  the  time  of  arrival  at 
the  second  shelter.  Thus,  the  time  interval  during  which  residents  are  shel¬ 
tered  as  a  result  of  the  third  visit  to  shelters  is  .  120  +  .  095  =  .  215  hours  and 


lA  ssuming  no  bottleneck  at  the  shelter  entrance. 
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.  143  +  .095  =  .238  hours  respectively.  This  procedure  is  continued  until  the  time 
intervals  have  been  established  for  all  residents  seeking  shelter  within  the 
tract.  The  results  of  applying  the  resident  sheltering  model  to  Tract  14  are 
contained  in  Table  2,  which  follows. 


Table  2 

Residents  Sheltered  as  a  Function  of  Time 
for  Census  Tract  14  of  the  City  of  Stamford,  Connecticut 


Residents  Sheltered 

Time  from  Departing  Homes 
(hrs. ) 

1300 

.  025 

2083 

.  032 

2315 

.  034 

2497 

.  036 

2935 

.  040 

3020 

.  041 

3444 

.  047 

3474 

.  048 

3474 

.  120 

3823 

.143 

3823 

.  215 

4063 

.  238 

4063 

.  310 

4192 

.333 

4192 

.  405 

4290 

.  428 

4290 

.  500 

4388 

.  528 

4388 

.  595 

4486 

.  618 

4486 

.  690 

4584 

.713 

4584 

.  785 

4682 

.  808 

4682 

.  880 

4781 

.  903 
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Model  to  predict  the  number  of  nontract  residents  sheltered 
within  a  tract  having  an  excess  of  shelter  spaces,  as  a  func¬ 
tion  of  elapsed  time  from  when  nontract  residents  arrive. 


Appendices  3  and  5  contain  two  models  which  represent  respectively: 

.  The  process  by  which  people  who  cannot  find  shelter  within 
their  own  tract  are  assigned  to  and  arrive  at  tracts  having 
shelters  with  excess  capacity.  The  end-product  of  this 
model  is  the  arrival  of  people  at  the  boundaries  of  tracts 
having  excess  capacity  as  a  function  of  time. 

.  The  process  by  which  people  who  are  located  within  a  census 
tract  enter  shelters  within  that  same  tract.  The  end-product 
of  thifr  model  is  the  number  of  people  sheltered  as  a  function 
of  time. 

The  model  described  here  picks  up  where  the  above  two  models  leave  off. 

It  predicts  the  number  of  nontract  residents  sheltered  within  a  tract  having 
an  excess  of  shelter  spaces,  as  a  fimction  of  time.  An  illustration  will  be 
used  as  an  aid  in  developing  the  model.  The  illustration  is  for  census  tract 
#12  in  the  City  of  Stamford,  Connecticut. 

I.  REQUIRED  DATA 


The  following  data  are  a  prerequisite  for  use  of  the  model. 

1.  The  number  of  shelters  in  the  tract  being  considered,  having  possible 
protection  factors  of  100  or  better,  and  the  remaining  capacity  of  each  shelter 
after  residents  of  that  tract  have  taken  shelter.  Table  1  contains  this  informa¬ 
tion  for  tract  #12,  City  of  Stamford,  Connecticut.  The  number  of  shelters  to 
be  considered  plus  the  "Total  Shelter  Spaces  Before  Tract  Residents  Take 
Shelter"  is  taken  from  the  National  Fallout  Shelter  Survey,  Phase  2,  under  the 
classification  "Total  Possible  Shelter  Spaces,  with  protection  factor  100  or 
better.  "  "Shelter  Capacity  After  3,  324  Tract  #12  Residents  Take  Shelter"  is 
determined  using  the  model  described  in  Appendix  4. 
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Table  1 


Shelter  Capacity  in  Tract  #12 


Facility  Number 

Total  Shelter  Spaces 
Before  Tract  Residents 
Take  Shelter 

Shelter  Capacity  After 

3,  324  Tract  #12 
Residents  Take  Shelter 

1 

353 

0 

2 

695 

0 

3 

2,625 

1,487 

4 

3,043 

1,905 

2.  The  area  of  shelter  tract:  For  tract  #12  the  area  is  1.055  square 
miles.  This  information  can  be  obtained  in  a  number  of  ways.  In  this  case 
a  planometer  was  used  on  a  1/9600  scale  map  of  Stamford. 

3.  The  number  of  people  from  other  tracts  seeking  shel^r  in  the  tract 
being  considered,  and  their  time  of  arrival  and  rate  of  arrival  at  the  tract: 
For  tract  #12  this  information  is  listed  in  Table  2. 

Table  2 

Number  and  Arrival  Rate  of  Nontract  #12 
Residents  Arriving  in  Tract  #12 


Number  of  Nontract 

#12  Residents  Arriv-  Time  of  Arrival  Time  Period  Rate  of  Arrival 


ing  in  Tract  #12 

(minu\ues) 

(minutes) 

(People /Minute) 

1355 

0  -  22.  8 

22.8 

59.43 

0 

22.8  *  41.4 

18.6 

0 

2037 

41.4  -  58.  1 

16.  7 

121.63 

The  data  in  this  table  are  determined  by  using  the  model  described  in 
Appendix  3. 
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IL  THE  USE  OF  ASSUMPTIONS  IN  THE  MODEL 


This  model  is  based  on  a  number  of  assumptions.  These  assumptions 
are  listed  below  along  with  appropriate  illustrations  of  their  use  in  pre* 
dieting  the  (Sheltering  of  a  nonresident  population  in  tract  #12  shelters. 

1.  The  shelters  within  a  tract  are  assumed  to  be  uniformly  distri¬ 
buted  within  that  tract  and  to  service  equal  areas  such  that  the  sum  of 
theee  areas  equals  the  area  of  the  entire  tract.  Under  these  conditions 
the  distance  between  shelters  can  be  represented  by  equation  (1): 


where:  d  =  the  distance  between  adjacent  shelters  in  miles. 

A  =  the  area  of  the  tract  in  square  miles. 

N  =  the  number  of  shelters  having  possible 
protection  factors  of  100  or  better* 

In  this  case.  A  =  1.  1  square  miles  (page  6o2)»  and  N  -  4  shelters  (page  6>2), 
therefore,  for  tract  #12,  d  -  .  52  miles. 

2.  It  is  assumed  that  each  shelter  within  a  tract  has  an  equal  likeli¬ 
hood  of  being  the  first  shelter  a  new  arrival  finds,  and  that  if  this  shelter 
is  full,  safety  is  sought  in  the  next  nearest  shelter.  All  movements  are 
straight-line  distances  from  one  shelter  to  another.  This  process  con¬ 
tinues  until  a  persou  arrives  at  a  shelter  having  excess  shelter  capacity. 

3.  All  movement  toward  shelters  within  a  tract  is  assn  ned  to  be  on 
foot,  at  the  rate  of  3  miles  per  hour  or  .05  miles  per  minute. 

4.  People  arrive  at  a  tract,  in  search  of  shelter,  either  on  foot  or 
by  automobile.  If  they  arrive  on  foot,  measurement  of  the  distance  they 
travel  to  their  "first"  shelter  begins  as  they  cross  into  the  tract.  If  they 
arrive  by  car,  the  distance  to  their  "first"  shelter  is  measured  from 
wherever  they  park  their  car  to  the  nearest  shelter.  In  both  cases  this 
initial  distance  to  a  "first"  shelter  can  be  represented  by  equation  (2): 

R  *  .38id  (2) 
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where:  R  =  mean  distance  to  the  first  shelter  in  miles. 

d  =  distance  between  adjacent  shelters  in  miles 
(equation  (1)). 

In  this  case,  d  =  .  52  miles;  therefore,  R  for  tract  #12  =.  20  miles. 

5.  All  shelters  with  excess  capacity  are  assumed  to  fill  up  at  the  same 
rate.  On  the  basis  of  this  assumption  it  is  possible  to  calculate  the  incre¬ 
mental  number  of  persons  who  must  take  shelter  before  each  additional 
shelter  is  filled.  For  the  example  used  here  this  situation  is  described 
in  Table  3. 


Table  3 

Number  of  Nontract  Residents  Required  to  Fill 
the  Third  and  Fourth  Shelters  in  Tract  #12 


Facility 

No. 

Shelter  Capacity 
After  3,  324  Tracks 
#12  Residents 

Take  Shelter 

Number  of  Nontract 
Residents  who  Must 
Take  Shelter  Before 
One  Additional  Shelter 
is  FUled 

Number  of  Remain¬ 
ing  Shelters  Having 
Excess  Capacity 

1 

0 

0 

Two 

2 

0 

0 

Two 

3 

1,487 

2.974 

One 

4 

1,905 

418 

None 

6.  It  is  assumed  that  a  person  seeking  shelter  will  visit  any  particular 
shelter  only  once;  i.  e. ,  having  once  visited  a  full  shelter,  he  will  not  return 
to  it.  The  consequence  of  this  assumption  is  that  the  universe  of  shelters  a 
person  might  visit  grows  smaller  each  time  he  finds  a  full  shelter.  For 
example,  suppose  a  person  seeks  shelter  in  a  tract  having  6  shelters  under 
the  conditions  where  there  are  4  full  shelters  and  4  shelters  with  excess 
capacity.  He  has  a  50%  chance  of  finding  a  shelter  with  excess  capacity 
\xpon  arrival  at  the  first  shelter.  If  this  first  shelter  happened  to  be  full, 
he  then  has  3  full  shelters  and  4  shelters  with  excess  capacity  to  choose 
from.  Therefore,  his  chances  of  finding  space  in  a  shelter  on  his  second 
try  are  4  out  of  7  or  57%. 
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7.  During  each  period  where  the  number  of  shelters  having  excess 
capacity  remains  constant*  it  is  assumed  that  the  mean  distance  a  new 
arrival  must  travel  from  the  tract  border  to  a  shelter  having  excess  capa¬ 
city  also  remains  constant.  This  mean  distance  is  represented  by  equa¬ 
tion  (3): 

D  =  R  +  ^(d)+(^)(|)(d)  +  (^)(~|)(^)(d)....  (3) 
where:  n  =  number  of  shelters  with  no  remaining  capacity. 

N  -  number  of  shelters  having  protection  factors 
of  100  or  better. 

R  =  mean  distance  from  the  tract  border  to  the  first  shelter 
in  miles  (equation  (2)*  page  h.3). 

d  =  distance  between  adjacent  shelters  in  miles  (equation 
(1),  page  6-3). 

O  =  mean  distance  from  tract  border  to  a  shelter  having 
excess  capacity,  in  miles. 

In  the  example  used  here  it  is  necessary  to  calculate  two  values  of 
D.  The  first  applies  to  that  group  of  people  seeking  shelter  under  the  initicl 
conditions  where  two  of  the  four  possible  shelters  have  excess  capacity.  In 
this  case: 

n  =  two  shelters  have  no  remaining  capacity. 

N  =  four  shelters  having  possible  protection  factors  of  100  or  bette 

R  =  .  20  miles. 

d  =  .  52  miles, 
and*  therefore* 

Dj  =  .55  miles. 

Similarly,  when  the  third  shelter  fills  up*  n  now  equals  three,  and  i' 
calculated  as  l.OO  miles. 


8.  Each  time  a  shelter  is  iilled»  it  is  assumed  that  no  additional 
persons  enter  shelter  for  a  period  of  time  represented  by  equation  (4): 

D-D 

^ _ 1  (4) 

r 

time  in  minutes  during  which  no  one  enters  shelter, 

mean  distance  in  miles  from  tract  border  to  an  empty 
shelter  immediately  following  the  filling  up  of  the  last 
shelter  {equation  (3),  page  f^-5), 

mean  distance  in  miles  from  tract  border  to  an  empty 
shelter  immediately  previous  to  the  last  shelter  filling 
up  (equation  (3),  page  6-5). 

rate  at  which  people  walk  to  shelter  within  the  tract 
(required  data  page  6-3). 

In  the  case  of  tract  #12,  this  situation  occurs  only  once,  at  the 
point  in  time  when  the  third  shelter  fills  up.  Under  these  conditions, 

Dy  =  Dj  -  .55  miles.  D^  =  D2  =  1.00  miles,  and  r  =  .  05  miles 
per  minute.  Consequently,  t  =  9.  00  minutes, 

IIL  APPLICATION  OF  THE  MODEL 


where :  t  = 

D  = 

X 


r 


Using  the  required  data  and  assumptions  presented  in  Sections  I  and 
II,  the  rate  at  which  shelter  spaces  are  filled  by  persons  arriving  from 
other  tracts  can  be  calculated.  The  step  by  step  process  by  which  these 
calculations  are  made  for  tract  #12,  City  of  Stamford,  Connecticut,  is 
presented  in  Table  4,  Combining  the  information  in  columns  "C"  (time) 
and  "E''  (total  people  in  shelter)  of  Table  4,  it  is  possible  to  create  a 
graphic  representation  of  nontract  residents  sheltered  as  a  function  of 
time.  Such  a  graph  for  tract  #12  is  presented  on  page  6-9,  Similar 
graphs  are  constructed  for  each  census  tract  having  a  shelter  capacity 
greater  than  the  resident  population.  For  the  City  of  Stamford,  Connec¬ 
ticut,  this  is  the  case  for  tracts  1,  11,  12,  16,  17,  18,  22,  and  23,  By 
combining  graphs  of  this  type  with  similar  graphs  showing  entrance  into 
shelter  as  a  function  of  time  for  tract  residents,  it  is  possible  to  develop 
a  composite  graph  for  an  entire  area  or  city.  Such  a  graph  for  the  City 
of  Stamford,  Connecticut,  is  presented  on  page  20  in  Volume  I, 
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Step  by  step  process  by  which  shelters  in  Tract  #12,  City  of  Stamford,  Conn.  , 
not  filled  by  Tract  #12  residents  are  filled  by  people  arriving  from  other  tracts. 
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Figure  1,  Residents  of  Other  Stamford  Tracts  Sheltered  in  Tract  #12  versus  Flapsed 
Time  from  Departing  Homes,  Assuming  40  mph  Driving  5peed,  3  mph 
Walking  Speed,  One  Lane  Traffic,  and  No  Individual  Shelter  Assignment. 
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APPENDIX  7 


Mathematical  Description  of  Fire  Probabilities 


1.  The  State  Eqrations 

It  is  postulated  that  initially  the  depopulated  region  contains  N  "sources" 
of  potential  fires:  typical  examples  are  unextinguished  cigarettes,  operating 
electrical  appliances,  automatic  heating  systems,  and  structures  inadequately 
protected  from  lightning.  During  any  time  interval  [  0,  t  j  one  of  three  events 
associated  with  such  a  source  may  occur. 

(a)  The  source  may  lose  its  potential  without  initiating  a  self- 
sustaining  fire:  the  cigarette  may  burn  out  harmlessly,  the 
fuse  in  the  circuit  feeding  the  appliance  may  blow,  etc. 

(b)  The  potential  may  be  converted  into  actuality:  the  source 
may  initiate  a  self-sustaining  fire  in  its  surroundings. 

(c)  The  potential  may  continue  unmodified:  the  heating  system 
may  operate  and  remain  a  potential  hazard,  the  danger  from 
electrical  storms  continues  unmodified,  etc. 

These  events  occur  at  random  points  in  time.  As  stated  in  the  main  text,  the 
assumption  used  in  this  development  is  that  the  probability  that  such  an  event 
occurs  in  an  interval  [t,  t  +  h],  given  that  nothing  which  makes  the  event  im¬ 
possible,  occurred  during  [0,t],  is  proportional  to  the  duration  of  the  interval, 
h,  with  some  undefined  dependence  upon  terms  of  order  h^  or  larger.  The 
formal  application  of  this  assumption  to  the  three  cases  listed  is  as  follows. 

(a)  Given  that  neither  an  event  of  type  (a)  nor  an  event  of  type  (b) 
occurred  during  [0,  t],  the  probability  that  the  event  of  type  (?,} 
occurs  in  [t,  t  +  h]  is  ah  +  o(h).  Here,  "a"  is  a  proportionality 
constant,  and  o(h)  represents  terms  of  order  h^  or  higher. 

(b)  Given  that  neither  an  event  of  type  (a)  nor  an  event  of  type  (b)  occurred 
in  [0,t],  the  probability  that  an  event  of  type  (b)  will  occur  in  [t,  t+h]  is 
bh+o(h),  where  "b"  is  a  constant  of  proportionality. 

(c)  It  follows  that  under  the  same  given  conditions  the  probability 
that  neither  type  of  event  occurs  in  [t,t  +  h]  is  1  -  (a  +  b)h  +  o(h). 
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It  is  convenient  to  describe  any  possible  combination  of  events  in  the 
region  by  defining  "states  of  the  region"  as  follows. 


state  0:  no  events  of  either  kind  for  any  of  the  N  sources 

state  1:  one  event  of  type  (a),  none  of  type  (b) 


state  n:  n  events  of  type  (a),  none  of  type  (b) 


state  N:  N  events  of  type  (a),  none  of  type  (b) 
state  Z:  one  event  of  type  (b),  number  of  prior  events  of 

type  (a)  not  specified. 


The  region  has  two  terminal  states.  If  it  is  in  state  N  no  further  events 
can  occur,  and  if  it  is  in  state  Z  (on  fire)  additional  events  of  type  (a)  are 
no  longer  significant. 


The  procedure  used  by  Feller*  is  as  follows.  Let  Pn(t)  represent  the 
probability  that  the  region  is  in  state  n,  n  ^  z,  at  time  t.  The  probability 
that  the  region  is  in  state  n  at  t  +  h  is  the  probability  of  the  following  com¬ 
pound  event: 


(a)  n  events  of  type  (a)  occurred  during  [0,t]  and  none  during 
[t,t+h];  the  probability  of  this  event  is  Pn(t)  [l  -  (a  +  b)h +o(h)]^"^ 
since  there  are  (N-n)  sources  available; 

(b)  (n  -  1 )  events  of  type  (a)  occurred  in  [0,t],  one  event  of  type  (a) 
occurred  in  [t,t  +h]  and  none  of  type  (b)  in  [t,t  +h]:  the  probability 
of  this  event  is  Pn  -  i(t)[ah +o(h)]. 

Accordingly  the  probability  of  the  compound  event  is 

(1)  Pn(t+h)  >  Pn(t)[l  -  (N-n)ch+o(h)]  +  Pn. i{t)[ah +o(h)] 

where  c  x  a  +  b  has  been  introduced  for  convenience. 


*Feller,  William,  An  Introduction  to  Probability  Theory  and  its  Applica¬ 
tions,  Vol.  I,  John  Wiley  and  Sons,  New  York;  1950.  In  particular,  see 
Chapter  17. 
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Rearranging  terms 


o(h) 

h 


P  (t+h)-P  (t) 

(2)  —2 - - - 2_  =  _(N-n)cP  (t)  +  aP  (t)  + 

n  n  n- 1 

In  the  limit  as  h  approaches  zero,  the  term  on  the  left  becomes  the  deriva¬ 
tive  of  Pjj  and  the  last  term  on  the  right  approaches  zero.  Thus 

(3)  P'(t)  +  (N-n)cP  (t)  =  aP  ,(t)  ;  n  =  1,  2,  .  .  . ,  (N- 1) 

n  n  n- 1 

where  the  prime  indicates  differentiation  with  respect  to  t.  Differential 
equations  for  PQ(t)  and  can  be  derived  by  an  analogous  process 

modified  to  take  account  of  the  fact  that  there  are  no  transitions  to  state 
zero  and  no  transitions  from  state  N.  The  resultant  equations  are 


(4)  Pyt)+  NcPQ(t)=  0 
and 

(5)  PJj(t)  =  aP^_j(t) 

(A  differential  equation  for  p2(t)  can  also  be  drived  in  this  way,  but 

a  simpler  procedure  is  to  compute  this  from  n?(fn^^^  "  The  initial 
conditions  to  be  satisfied  by  these  equations  are  derived  from  the 
assumption  that  no  events  occur  before  t  =  0:  these  are  P  (0)  =  1,  and 
Pjj{0)  =  0.  It  is  shown  in  Supplement  A  to  thisSection  that  the  solutions  to 
these  equations  are 

n 

(6)  P(t)=4j-  exp  C-(N-n)ct] 

n  ni 

where 

(7)  r  =  [1  -  exp{-ct)] 

For  later  reference  it  is  noted  that 

(8)  r '  =  a  exp(  -  ct) 
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2.  The  Physical  Significance  of  the  Proportionality  Constants 


In  the  preceding  discussion,  the  constants  a,  b,  and  c  =  a  f  b,  were 
introduced  as  constants  of  proportionality.  Considered  dimensionally, 
these  constants  represent  rates;  that  is  a  number  per  unit  of  time.  Refer¬ 
ring  back  to  the  definition  of  "a*,*  it  can  be  seen  that  if  this  is  increased, 
the  probability  of  a  non- fire  producing  event  during  a  given  interval  is 
increased.  It  follows  that  "a"  represents  the  rate  at  which  such  events 
Occur.  Similarly,  •'b"  represents  the  rate  at  which  fire-producing  events 
occur.  A  more  precise  characterization  of  the  significance  of  these  con¬ 
stants  can  be  obtained  as  follows. 


Consider  a  region  containing  only  one  source;  i.  e. ,  N  =  1.  Let  the 
time  scale  be  divided  into  intervals  of  length  h  so  that  tj  =  ih.  The  life 
(or  duration)  of  state  0  can  be  assigned  the  value  tj  if  no  events  occur 
during  [O,  tj]  and  an  event  of  type  a  or  type  b  occurs  during  [tj,  tj  +  h]. 
Since  the  probability  of  this  event  is  ch  the  expected  value  of  the  state  0 
lifetime  is 

OD 

(9)  T  =y  t.P  (t.)ch 

o  L  1  o  1 
o 

In  the  limit,  as  h  approaches  zero. 


(10) 


T 

o 


00 

n 


00 

p 


tP  (t)  dt  =  c 


t  exp(-ct)dt  =  (1/c) 


0 


0 


If  an  ensemble  of  regions  with  N  =  1  is  considered,  Tq  can  also  be  called 
the  mean  accident-free  time  and  the  ergodic  implication  is  that  (1/c)  can 
also  be  called  the  mean  time  between  accidents  if  the  state  is  restored  to 
zero  after  each  accident.  It  is  also  consistent  to  interpret  c  as  the  average 
accident  rate. 


The  same  rationale  leads  to  the  conclusion  that  the  expected  value  of 
the  time  required  to  reach  the  state  P^,  after  which  no  accident  is  possible, 
is  given  by 


(11) 


T  =  a 


tP  (t)dt 
o 


(a/c^) 


and  the  expected  value  of  the  time  required  for  a  fire  to  start  (i.  e. ,  to 
reach  state  Z)  is  given  by 
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(12) 


r  2 

T  =  b  tP  (t)dt  =  (b/c  ) 

Z  J  o 

0 

In  other  words,  (c^/b)  is  the  average  rate  at  which  the  accident  potential 
is  dissipated  harmlessly,  and  (c^/a)  is  the  average  rate  at  which  fires 
start. 

The  rationale  can  be  extended  to  the  general  case.  State  Z  can  be 
reached  by  a  transition  from  any  state  except  n=  N.  Since  the  probability 
of  a  tran^ion  from  any  state  n/N  to  state  Z  is  (N  —  n)bh,  the  expected 
value  of  the  time  before  fire  occurs  is 

N-1  f 

(13)  T_  =  b  r  (N-n)  tP  (t)  dt 

n=0  J 

where  the  order  of  integration  and  summation  have  been  interchanged. 
This  quantity  is  evaluated  in  Supplement  B,  giving 

(14)  T  =  b/Nc^ 


Alternately,  Nc  /b  is  the  average  rate  at  which  fires  start  in  an  ensemble 
of  regions  with  N  sources  in  each  region. 


3.  The  Cluster 


The  preceding  discussion  has  described  events  in  a  region  containing 
N  potential  sources  of  spontaneous  ignition,  and  which  burns  completely 
if  one  of  these  sources  starts  a  self-sustaining  fire.  It  is  convenient  to 
refer  to  such  a  region  as  a  cluster:  this  name  suggests  that  the  structures 
and  ignition  points  are  distributed  in  a  way  which  makes  fire  spread  through 
the  region  inevitable  if  no  control  is  attempted.  A  more  complete  physical 
characteristic  of  a  cluster  must  await  further  developments  by  the  experts 
studying  fire  spread  mechanisms.  It  is  clear  that  a  cluster  may  be  an 
entire  downtown  region  or  a  single  house  in  the  suburbs.  The  utility  of  this 
concept  will  become  apparent  later  in  the  discussion  when  it  is  pointed  out 
that  fire  spread  can  be  characterized  by  considering  the  probability  that  a 
burning  cluster  ignites  a  contiguous  one.  The  statistical  description  of  such 
events  can  be  developed  despite  the  lack  of  a  clear  physical  definition. 
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4.  The  State  Equations  for  a  Cluster 


The  condition  of  a  cluster  can  be  characterized  by  defining  three 
states  as  follows: 

State  0  -  no  fire  has  occurred 
State  1  -  cluster  burning 
State  2  -  fire  burned  out 

It  is  recognized  that  the  words  "burning"  and  "out"  used  in  defining  States 
1  and  2  are  not  precise.  In  a  more  exact  analysis  than  is  of  interest  here, 
it  would  certainly  be  appropriate  to  introduce  the  states  "violent  burning" 
and  "residual  burning"  discussed  by  Chandler,  et  al.  The  discussion  is 
materially  simplified  by  ignorning  this  distinction  for  the  moment. 

Using  this  definition  of  the  states,  the  events  of  interest  can  be 
characterized  as  transitions  from  one  state  to  another.  The  transition 

0 - ►!  represents  the  event:  the  cluster  catches  fire.  The  transition 

1 - ^2  represents  the  burning  out  of  the  fire.  Clearly,  there  are  no 

transitions  from  State  2. 

Using  these  assumptions  and  definitions  and  the  analytical  procedure 
described  in  Section  4,  the  following  results  are  obtained. 

(15)  Pg(t)=e‘*‘ 

describes  the  probability  that  the  cluster  has  not  caught  fire  by  time  t 
given  that  it  was  not  on  fire  at  time  0.  The  probability  that  the  cluster 
is  burning  at  any  time  t  is  given  by 

(16)  P  (t)  =-^  (e'“‘-e''’') 

1  b  —  a 

The  analogous  procedure  for  P^^t)  gives 

(17)  (be’^‘-ae'^^^ 

As  a  check  it  is  noted  that 

(18)  P^d)  f  Pj(t)  f  P^d)  -  1 


♦Chandler,  C.  C.  ,  Storey,  T.  G. .  and  Tangren.  C.  D.  Prediction  of  Fire 
Spread  Following  Nuclear  Explosions.  U.  S.  Forest  Service  Research 
Paper  PSW-^,  1963. 
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5.  Physical  Interpretation 


Using  the  established  definitions,  the  probability  that  the  transition 
from  0  to  -  Occurs  in  the  interval  [t,  t+  h]  is  ahPQ(t).  The  expected  value 
of  the  time  at  which  this  transition  takes  place  is 

n 

(19)  E(T  )  =  a^ 

i=  0 


Converting  this  to  an  integral  by  letting  n  approach  .nfinity  and  the  largest 
hj  approach  zero. 

00 

(20)  E(Tq)  =  a  fpQ(t)dt=  (1/a) 

0 

Thus  a  is  the  reciprocal  of  the  time  before  fire  starts  or  of  the  time  the 
cluster  is  in  State  0.  If  an  ensemble  of  independent  clusters  is  considered, 
this  can  also  be  interpreted  as  the  mean  (or  average)  lifetime  of  the  fire- 
free  State  0.  Alternately,  "a"  itself  can  be  interpreted  as  the  average  rate 
at  which  fires  start  in  the  ensemble  of  clusters.  IL  Is  clear  that  any  action 
taken  before  depopulation  (e.  g. ,  inspection  to  reduce  the  number  of  potential 
ignition  points  in  the  clusters)  will  decrease  a  or  increase  (1/a). 

An  analogous  analysis  yields  that  the  expected  value  of  the  time  at 
which  the  transition  from  State  1  to  State  2  occurs  is  giv«n  by 


(21) 


n 


E(T,)  =  bJ 


tPj(t)dt  = 


a  -t-  b 
ub" 


0 


The  difference  between  these  two  quantities  is  the  expected  value  of  the  time 
between  the  two  transitions,  or  the  expected  value  of  the  lifetime  of  the 
burning  State  1.  This  is 


(22)  a  f  b  ^  ±  ± 

ab  a  b 

Thus  b  is  the  reciprocal  of  the  average  duration  of  the  fire.  It  can  also  be 
interpreted  as  a  burn-out  rate  but  this  doc*  not  appear  to  be  very  helpful 
There  is  little  or  nothmg  that  can  be  done  before  depopulation  which  may 
affect  b  materially,  b  would  appear  to  be  a  (unction  of  fuel  density  and  the 
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size  of  the  cluster  with  some  possible  statistical  variations  induced  by  rain 
or  its  absence,  wind  conditions,  etc.  Because  of  this,  a  convenient  way  to 
normalize  the  description  is  to  let  (1/b)  equal  one  arbitrary  unit  of  time. 
Some  typical  results  when  b  =  1  are  shown  in  Figure  2  of  the  text. 


6.  A  Rough  Estimate  of  the  Hazard 

Finally,  it  is  interesting  to  speculate  about  the  possible  values  of  a 
and  b.  If  it  is  recalled  that  a  can  be  interpreted  as  the  mean  rate  at  which 
fires  occur,  an  estimate  of  its  value  for  normally  occupied  areas  can  be 
obtained  from  statistics  on  the  number  of  fires  which  actually  occur.  As 
a  typical  example,  the  number  of  actual  fires  in  Los  Angeles  in  the  year 
starting  1  July  1958  is  given*  as  13,  933.  For  this  case,  a  =  1.  6  per  hour, 
or  equivalently  ( l/a)  the  mean  time  between  fires  is  0.  625  hours.  We  have 
not  found  any  useful  guide  for  extrapolation  from  this  number  to  a  figure 
which  might  apply  during  depopulation.  Even  an  optimistic  doubling  of  this 
mean  time  however  leads  to  the  conclusion  that  this  is  a  serious  hazard 
which  may  render  any  other  measures  for  protecting  the  population  relatively 
unimportant,  especially  if  the  fallout  shelters  themselves  are  threatened  by 
the  spreading  fires. 

7.  An  Initial  Discussion  of  Fire  Spread 

The  analysis  developed  in  Section  4  can  be  extended  to  include  consid¬ 
eration  of  fire  spread  by  introducing  a  modified  definition  of  state,  and  using 
the  concept  of  a  cluster.  The  essential  characteristic  of  a  cluster  is  that 
there  is  a  finite  probability  that  it  will  burn  out  without  spreading  to  neighbor¬ 
ing  areas.  Conceptually,  any  region  can  be  subdivided  into  M  such  clusters. 
The  state  of  the  region  can  then  be  represented  as  a  M-tuple  (ij,  i^,  .  .  .  ,  ij^), 
in  which  each  element  represents  the  state  of  the  corresponding  cluster.  If 
the  cluster  states  defined  in  Section  5  are  used,  i  =  0  indicates  no  ignition, 
i  =  1  indicates  active  burning,  and  i  =  2  indicates  burned  out,  and  there  are 
then  3^  possible  states.  (More  detailed  definitions  of  cluster  states  based 
on  the  N  sources  in  each  cluster  can  also  be  used.  )  Both  initial  ignition 
and  spread  between  clusters  can  then  be  described  analytically  as  transitions 
from  one  of  these  states  to  another. 


^'A  Study  of  Fire  Problems,  Publication  949,  National  Academy  of  Sciences, 
National  Research  Council,  Washington,  D.  C. ,  1961.  (This  is  a  report  of 
a  study  held  at  Woods  Hole,  Massachusetts,  under  the  auspices  of  the 
Committee  on  Fire  Research  of  the  National  Academy  of  Sciences.  ) 
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As  a  comprehensible  example  which  can  easily  be  generalized,  suppose 
that  M-  2;  i.  e. ,  thatthe  region  is  divided  into  two  clusters.  Then  the  state 
(0,  0)  designates  that  no  fires  have  started  in  either  cluster.  The  transition 
(0,  0)  to  (0,  1)  indicates  that  a  fire  has  started  in  the  second  cluster  which  is 
now  burning.  The  probability  of  such  an  initial  transition  is  that  discussed 
for  transitions  from  0  to  1  in  a  single  cluster;  i.  e. ,  a2h  +  o(h)  where  the 
subscript  identifies  the  cluster.  A  subsequent  transition  may  be  from  (0,  1) 
to  (1,  1).  This  can  occur  in  one  of  two  ways:  there  may  be  spontaneous 
ignition  in  cluster  1,  or  the  fire  may  spread  from  cluster  2  to  cluster  1. 

This  second  possibility  can  be  described  by  introducing  a  probability  of  the 
form  gijh  +  o{h),  where  g  is  completely  analogous  to  the  parameters  a  and 
b  used  in  the  preceding  discussion.  The  probability  of  the  transition  from 
(0,  1)  to  (1,  1)  in  [t,  t+  h]  is  then  (aj  +  g2l)h  0(h)*  Other  possible  transi¬ 
tions  can  be  analyzed  in  the  same  way.  When  higher  order  terms  are  ignored, 
the  possible  transitions  are  those  indicated  by  the  arrows  in  Figure  2. 

When  the  spread  of  fire  is  considered  in  this  way,  it  is  clear  that  the 
methods  described  in  Section  4  are  directly  applicable  to  the  derivation  of 
simultaneous  linear  differential  equations  describing  the  probability  Pij(t) 
that  the  region  is  in  any  of  the  nine  possible  states  at  any  time  t.  From 
this  description,  it  is  also  clear  that  if  the  description  is  quantized,  as 
would  normally  be  done  if  a  digital  computer  is  to  be  used  for  the  calcula¬ 
tions,  the  process  is  a  simple  Markov  process;  namely,  a  random  walk  with 
one  absorbing  state,  (2,  2).  The  matrix  of  transition  probabilities  for  this 
particular  example  is  of  the  form  shown  in  Figure  3. 


w(0>  2) 

(0.  1) 

(0, 


Fig  ure  2.  Possible  transitions  for  region  with  two  clusters. 


(1,  1) 


(1.  2) 


"^2.  2) 
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Figure  3.  Matrix  of  transition  probabilities  for 
region  with  two  clusters. 
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SUPPLEMENT  A  TO  APPENDIX  7 


Solution  to  State  Equations 

Tb  demonstrate  that  equation  (16)  provides  solutions  to  equations  (3),  (4), 
and  (5),  it  is  sufficient  to  substitute  the  indicated  quantities  into  the  equations 
and  check  that  the  result  is  an  identity.  The  asserted  solutions  are 


(6) 


r  -(N  -  n)  ct 
P  =  — —  e 
n  nl 


where 

(7) 


r  = 


f(‘-‘ 


and 

(8) 


dr 

dt 


=  ae 


-ct 


For  n  =  o,  equation  (6)  gives 


(A-1) 


P  =  e 
o 


-Net 


differentiating 


(A-2) 


P'  =  Nee 
o 


-Net 


these  clearly  satisfy  equation  (4). 
For  n=N,  equation  (6)  gives 


(A-3) 


P  =  r 
N 


N 


Nl 


differentiating 

•’k  =  (S. 


A7-1 


For  n  =  N  -  1,  equation  ( 6)  gives 


{A-5) 


N-1 

■p  — 

^^N-1  “  {N-l)l 


-ct 

e 


Since  these  two  are  equal,  equation  (5)  is  satisfied. 

To  check  that  equation  (3)  is  satisfied,  equation  (6)  is  differentiated, 
giving 

/*  /V  X..I  i  ~(N-n)ct  -  -ct  , 

(A-6)  ~  "nT  ^  ^  -(N-n)cr] 

Also,  from  equation  (6) 

(A-7)  (N-„)cP_^=  ^ 


Adding  these  as  indicated  by  equation  (3)  yields  the  quantity 

(A-8)  - -  n-1  -(N-n+ l)ct 

^  (n-l)l  ® 

Since  this  is  a  from  equation  (6),  equation  (3)  is  satisfied  for 


n  =  1,  2,  .  ...  N-1 

Finally,  it  can  be  checked  by  inspection  that  the  initial  conditions  are 
satisfied. 
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SUPPLEMENT  B  TO  APPENDIX  7 


Mean  Time  to  Ignition  in  Regions  with  N  Sources 


T 

z 


is  given  by  equation  (13);  i.  e. , 


(13) 


N-1 


00 

o 


E  (N-n) 


n=0 


0 


tP  (t)dt 
n 


It  will  be  shown  that  it  is  only  necessary  to  develop  the  first  few  terms  of 
this  summation.  Define 


(B-1)  l(n;N)  =  (N-n)  c 

In  particular,  for  n  =  0 

(B-2)  1(0;N)  =  Nc 


tP  dt 
a 


tP  dt 
o 


0 


Using  equation  (3)  for  P 


(B-3) 


1(0:N)  =  - 


tP'dt 


Integrating  by  parts 

00 

00  r 

(B-4)  1(0;N)  =  [“tP^]Q 

C 

The  first  term  vanishes  since  P  (®)  =  0.  Expressing  P  explicitly  by  using 
equation  (B-5)  for  n  -  0. 


(B-5) 

The  proiess 


1(0;N) 


00 


0 

is  repeated  for  n 


I 

Nc 

1. 
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(B-6) 


GO 


I(1;N)  =  (N-  l)c 


tPjdt 


0 


Using  equation  (4)  for  n=  1 
(B-7)  I(l;n)  =  a 


tP  dt  - 
o 


tP'  dt 
1 


The  first  term  is  {a/Nc)I(o;N)  and  is  evaluated  from  {B-5).  The  second 
term  is  integrated  by  parts  using  that  as  t  =  0,  Pj(t)  =  0,  exponentially 
dominating  t.  This  yields 


(B-8) 


1(1  :N) 


+  P  dt 

k  A 


Using  equation  (4)  with  n  =  1 


(B-9) 


(N-  l)c 


P  dt  =  a 
^  u 

0 


P  dt  - 
o 


P^dt 


The  second  term  vanishes  since  P.(0)  =  P  (®)  =  0.  The  first  term  was 
evaluated  in  obtaining  (B-5).  Thus^ 


00 


(B-10) 


o 


Pjdt  = 


0 


N(N-  l)c 


Substituting  this  into  (B-8) 


(B-11) 


I;i;N) 


( 2N  -  1  )a 

2.  Z 

N  (N-  l)c 


Repeating  the  process  for  n  Z 


(B  12) 


I(2:N) 


(■IN'^-bN  t  2)a^ 
Z)c^ 


F rom  {l^)  and  ( B  - 1 ) 


(B-i3) 


/.  c; 


1(0. N)  t  ftl.N)  f  I(^.N)  ^  . 
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Using  the  values  obtained  in  (B- 5),  (B-ll),  and(B-12) 


(B-14) 


b  (2N-l)a 

z  ^2  N{N-l)c 

Nc 


(3N^-  6N  +  2) 
N(N-  1}^(N-  2) 


+  .  .  . 


Since  a  <  c,  and  the  denominators  are  increasingly  higher  order  than  the 
numerators,  the  higher  order  terms  can  be  neglected. 


A  first  approximation,  therefore,  is 


